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A computer simulation of some power plants, transmission 
systems, and propellers has been developed. 

Theoretical models have been established mainly from 
published ».odel s for the following components. 

1. prime movers 

a. aircraft derivative gas turbines 

b. lov/ speed diesel engines 

c. medium speed die el engines 

2. transmissions 

a. direct drive (low-speed diesel) 

b. reduction gear with and wit. -out reversing 
clutch ar r an ge me n t 

3. propulsors 

a. conventional fixed pitch propellers 

b. controllable reversible pitch propellers 
The computer simulation utilized the Dynamic System 

Simulation Program (BYSYS) on an IbK 1130 computer. jYSYS 
v;as developed by the Department of Mechanical Engineering 
at M.I.T. Representative accelerations and crashback 
maneuvers were computed ^sing toe various propulsion systems. 
Several critical situations v/ere established involving 
the sequencing of power plant, transmission, and propulsion 
transients. 
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INTRODUCTION 



The rapid advances made in the field of automatic controls since 
World War II have found more and more applications in the area of ship 
propulsion. One of the principal driving forces behind this is the 
desire to operate plants safely and at or near their "optimal*’ operating 
point. The safe coordinated operation of some ship power plants con- 
sisting of high performance prime movers, sophisticated couplings, and 
propellers is no trivial matter. This is especially true during periods 
when fast ship response is desired, e.g., emergency stopping or acceler- 
ating , 

To be able to control a propulsion plant, one must be able to 
model it and know how it behaves under various disturbances or changes 
in inputs. The alternative is a costly experimental program for each 
system. The problem for the control engineer now becomes one of developin 
or choosing a model that suitably describes the characteristics of the 
particular power train components. These .component models may be either 
physical or functional. A physical model may accurately describe a 
component’s characteristics using, for example, Newton’ s Laws of Motion, 
thermodynamic relationships, conservation of mass relations, etc. These 
provide a valuable insight of the physical processes occurring in the 
operation of a component, and enable one to perceive what simplifications 
can be made to a model. Some components, on the other hand, may prove 
to be very non-linear and to describe their operational characteristics 
using physical laws may result in mathematical expressions which are not 
suitable for use in a simulation study. In this situation a simpler 
mathematical expression is developed which closely approximates the 
component’s characteristics over a reasonable range of inputs. 



INTRODUCTION (Confd) 

The philosophy behind the study described in this thesis has been 

to keep the models simple, but yet reasonable. To quote ?. J. Kiviat, 

who prepared a RAND memorandum on digital simulation modeling concepts, 

"the model should only be as detailed as is necessary to 
answer the question at hand; it should be so designed, 
however, so it can be expanded to include more detail 
without inordinate cost in those model areas which have 
a high probability of becoming subsequent subjects of 
concern. "(1) 

Once the component models are selected, drive train models can be 
synthesized and prepared for simulation. No discussion will be made 
concerning the relative merits of digital simulation compared to analog 
or analog-digital (hybrid) simulation. The simulation literature is 
replete with these discussions. Digital simulation was decided at the 
outset. 

A major part of this study has been concentrated on writing computer 
programs for prime movers, propellers, and transient operations. These 
programs were developed to investigate various control schemes for marine 
power plants. 

The gas turbines selected for simulation are the Pratt and Whitney 
FT4A-2 and the General Electric LM2500-A. A medium speed and a low speed 
diesel were also selected. The medium speed diesel selected is the 
Mirrless KV Major 12; the low speed diesel selected is the Burmeistor 
and Wain 7K98FF. The main reasons for selecting these particular engines 
were the availability of suitable data at reasonable power levels. 
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1. GAS TURBINE MODELS 

1.1 Introduction 

Gas turbines, as applied to ship propulsion, usually connote 
the inclusion of a power turbine or •'free'* turbine, which 
provides the driving torque. The power turbine is driven by 
the hot exhaust gases of the "gas generator." Thus when 
marine engineers refer to gas turbines, they usually include 
the combination of gas generator and power turbine in their 
meaning. 

There are, of course, exceptions to the statement con- 
cerning separate power turbines. The Koelin Class escoi’t 
frigates built for the German Federal Navy are driven by 
gas turbines which do not have a separate power turbine. 

The gas generator is coupled through a clutch to the drive 
train .(2) 

1.2 Aircraft Gas Turbine Models 

During the literature survey part of this study, several 
gas turbine models useful for aircraft propulsion simulation 
were found. In Reference (3), for example, Cottington pre- 
sents a simulation study of a two-spool turbojet engine, and 
also covers the design of a digital three-terra controller 
using digital simulation methods. In trying to determine an 
"optimum" set of gains, he treated the high pressure rotor 
system as a critically damped second order system. Values 
of these optimized gain settings for varying values of high 
pressure rotor speeds are presented. 

Mueller has developed a single spool gas turbine model 
from one-dimensional flow relationships. (4) He assumes that 
the only important dynamics are those of tho rotating compressor- 
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1. GAS TURBINE MODELS (Cont'd) 

1.2 Aircraft Gas Turbine Models (Cont'd) 

turbine spool. Mueller justifies this by noting that gas 
disturbances travel through the engine at the speed of sound, 
whereas a spool acceleration or deceleration takes some seconds. 
No results of his simulations are presented. 

In References (5) and (6), models are derived from a 
linear analysis of gas turbine dynamics. The model presented 
in the latter reference is highly detailed and is developed 
using non-dimensionalized parameters. To use this model, 
one must have access to the machinery performance maps. The 
prominent feature of this model is that it assumes stepwise 
linearity. The compressor performance map is utilized to 
obtain the coefficients in the linearized equation at each 
operating point. 

In References (7) and (8), Saravanamuttoo and Fawke dis- 
cuss the application of both analog and digital computer methods 
of simulating gas turbine dynamic performance. They approached 
the problem from the viewpoint of engineering thermodynamics, 
using the normal compressor and turbine characteristics. 

Figure (1-1) shows the results obtained by Saravanamuttoo and 
Fawke for a single-spool turbojet simulation, both digital 
and analog, along with the results of actual engine tests. 

It can be seen that the acceleration predicted by both simu- 
lations was much faster than the actual engine. An examination 
of the test results on the engine showed that there was a lag 
in the rise of turbine entry temperature following a fuel change. 
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1. GAS TURBINE MODELS (Conf d) 



1.2 Aircraft Gas Turbine Models (C.»r.t*d) 

Saravanamattoo and Fawke’s model, however, had assumed 
Immediate response, although they realized there would be 
some lag due to the dynamics of the fuel injection system. 

They decided to simulate this lag by postulating that after 
a sudden change in fuel flow, the combustion efficiency drops 
instantaneously and then recovers exponentially. This drop 
in efficiency was made proportional to the percentage increase 
in fuel flow. The constant of proportionality \ms determined 
such that the turbine entry temperature before and after the 
transient was unchanged from the steady state values. The 
time constant for the recovery of combustion efficiency was 
taken to be inversely proportional to the air flow through 
the combustor, i.e., 

Z = Z .Ha 

6 W 



Z - time constant for recovery of combustion efficiency 
- time constant at design speed 
air flow at design speed 
W - air flow 

To achieve flow compatability in their digital model, 
Saravanamattoo and Fawke added lags in the propagation of 
pressure signals through the engine. The lags were intro- 
duced by considering the engine components to be separated 
by specified volumes; these volumes produced the desired 
capacitive effect. 
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1. GAS TURBINE MODELS (Cont'd) 

1.3 Karins Gas Turbine Models 

During the early parts of the literature survey, a Russian 
title on analog and digital models of marine gas turbines 
was found in a gas turbine newsletter. This book, listed 
as Reference (9)» was written by Bichayev and published by 
the Shipbuilding Press of Leningrad in 19^9 ♦ Although it 
is still untranslated, the book appears to bo a major work 
in the field of marine gas turbine simulation, judging from 
the many analog computer diagrams, compressor maps, and 
equations. The cost of having it translated (over $900) 
was prohibitive. 

An unedited rough translation of a Russian work by Lurye 
discusses mathematical modeling of a marine gas turbine engine 
with a real fuel apparatus .(10) Lurye ‘ s model is presented 
in the form of linearized differential equations. The equations 
describing Lurye' s model are reproduced in Appendix A. 

Using Saravanamuttoo ' s technique, one can extend his 
model to a gas turbine having a free turbine. This is also 
presented in Appendix A. 

1 .4 Torque-Speed. Characteristics Model 

C. J. Rubis has developed a method based on the free 
turbine torque-speed characteristics to model a gas turbine's 
action. (11,12) The free turbine torque-speed characteristics 
can be derived from the steady-state power-speed characteristics 
which show lines of constant SFC. The corresponding curves 
for the Pratt and Whitney FT4A-2 and General Electric LK2500-A 
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1 • GAS TURBItJE MODELS (Ccnt«d) 

1*4 Torque-Speed Characteristics Model (Cont*d) 

are given in Appendix A. The torque-speed characteristics 
have been used here as the model selected for simulating the 
performance of the marine gas turbine. 

The unrefined model, used to describe the transient 
performance of the turbine operates in the following manner. 
Referring to Figure (1-2), assume that the initial steady-state 
operating point is at A operating at torque Tq and speed Nq* 
and assume that the engine is to accelerate to point B« In 
a time interval At, fuel is increased by A F. This increase 
in fuel produces an instantaneous change in torque A T such 
that the new torque is Tq + AT. This increase in torque will 
result in a change in power turbine rotation A N along the 
constant fuel rate line Fq + A F. The new starting point for 
the next time interval is at the intersection of the F 0 +AF 
constant fuel rate line and the new rotational speed Nq +A N* 
The dotted lines illustrate the paths taken during the time 
steps, while the heavy line is the torque path. 

1.5 An Improved Model 

The methodology of this model is not nextf. The torque- 
speed characteristics have been used in the past to model 
diesel engines. The idea of obtaining instantaneous torque 
for a change in fuel rate is much closer to the real thing in 
a diesel than it is in a gas turbine. Comparisons with engine 
tests have shown that the predicted values of torques obtained 
from the gas turbine characteristics are higher than the 
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1. GAS TURBItJE MODELS (Cont»d) 

1.5 An Improved Model (Cont*d) 

measured values. The differences in torque and the resulting 
differences in speed were measured by Rubis who conducted 
engine tests with the FT4A-2. He also used the torque-speed 
curves in a simulation between the same two power levels. 

Figure (1-3) shows how the simulation results in a higher 
acceleration because of the higher torques. This graph should 
be compared with figure ( 1 -/). 

It is evident from the comparison of figures (l-f) and (i-2>) 
that what is needed to improve this model is a correction 
factor, which could have a physical explanation similar to 
that offered by Saravan smut too and Fawke for their model; 
namely a reduction in combustion efficiency. The results of 
other transient tests conducted with the FT4A-2 were presented 
by Rubis in Reference (13). In this paper, Rubis presents 
correction factors to the torque-speed curve model for engine 
accelerations. This curve is reproduced herein as figure (1-4). 
One can observe that the curve for accelerations from idle can 
be expressed as a decaying exponential similar to the proposal 
by Saravansrauttoo and Fawke, viz . , 

Ke = .3 + .7 (l~e ) 

This expression, which creates a form of delay or lag in the 
transient torque produced by the model, was included in the 
simulations in this study. It should be clear that this lag 
is applied to the transient torque obtained from the model 
and not to the -whole torque. 
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1. GAS TURBINE MODELS (Cont»d) 

1*5 An Improved Model (Confd) 

The model appears to bo acceptable for ship propulsion 
studies because it is simple; it also permits loading other 
characteristics of the ship drive train into a medium sized 
computer without exceeding its core size. 
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2. D IESEL E NGINE 



2.1 Introduction 

The traditional use of steam power plants to propel 
large American built ships has created somewhat of a void 
in this country’s technical literature on diesel ship pro- 
pulsion. The successful use of diesel engines for ship 
propulsion by foreign shipbuilders, on the other hand, has 
been an established fact long before the first digital com- 
puter was ever built, and long before analog and digital 
simulations became established as accurate techniques for 
modeling and testing ship propulsion plants prior to their 
construction. .Consequently there has been very little demand 
from industry, either foreign or domestic, for diesel plant 
simulations. However , some published data were obtained 
and used in this study. Since the use of high speed diesels 
has been fairly restricted to small craft, the simulations 
have been limited to medium and low speed diesels. 

2.2 Presentation of Models 

Woodward presents in Reference (14) a simple, straight- 
forward treatment of the diesel engine as applied to ship 
propulsion. A widely used expression given in Reference (14) 
and other engine textbooks, is the formula for indicated 
horsepower (IHP): 



™ p = §io (HP) 

P = indicated mean effective pressure (lbs/in ) 

L = stroke (ft) 

O 

A = total piston area (in') for a multi-cylinder engine 
N = RPM (two-stroke engine), (four-stroke engine) 
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2. DIESEL ENGINE (Cent* d) 

2.2 Presentation of Models (Cont*d) 

This power is that produced in the engine cylinders, and is 
greater than that available at the output shaft because of 
friction and other internal loads (e.g., viater pumps, fuel pumps, 
turbochangers , etc.). What is available at the output shaft 
is termed brakehorsepower (BHP). The ratio of brakehorsepower 
(BHP) to indicated horsepower (IKP) is the mechanical efficiency 

(T^) of the engine, i.e., 

m = BH£ 

<m HIP 

or BHP = W PUN 
« n 

33$xT 

The product N [fl P is called the brake mean effective pressure 
(BMEP). The indicated mean effective pressure, P, can be 
directly measured and is proportional to the energy released 
by combustion in the cylinder. This energy is proportional 
to the amount of fuel injected in each power stroke. Most 
diesels utilize positive displacement pumps for fuel injection. 
The amount injected per power stroke for this type of pumping 
system should be the same. The stroke, L, and the total piston 
area, A, are, of course, fixed for a particular diesel. If 
one assumes constant mechanical efficiency, then it can be 
seen that for a particular fuel setting, BHP. is a linear 
function of N, i.e., 

BHP - KN 

K = N PIA for a particular fuel injection setting 

33000 
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2 * DIESEL ENGINE (Confd) 

2.2 Presentation of Models (Cont’d) 

Since rotational power is the product of torque, £, and speed, 
dc», it then follows that for a two stroke engine, 
co = 2 IT N (rad/min) 

Z<v = Z 2 ff N (ft- lb/ min) 

2^ = K^_ K as previously defined 

W 

i.e., torque is constant for particular fuel injection setting. 
The above is known as the ideal torque-rpra and power-rpm 
diesel relationships, and are graphically illustrated in 
Figure 2-1 and 2-2. 

The important assumption of this model is the one concern- 
ing constant mechanical efficiency. How close does this 
approach the actual case? In Reference (15), Pounder devotes 
an entire chapter to the Burraeister and Wain Engines and states 
a mechanical efficiency of 91$ in his discussion on B &W low 
speed diesels. (Note: this figure was confirmed by a represents 

tive of the Burmsister & Wain Corporation on a recent visit to 
MIT). 

It was indicated before that the indicated mean effective 
pressure is proportional to the fuel injected rate. An ex- 
tension of the ideal model then is to say that torque is not 
only constant for a given fuel setting, but also proportional 
to that fuel setting. In other words, torque is a linear 
function of fuel setting. Obviously, other extensions can 
be made such as relating torque to fuel rack position or to 
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2* DIESEL BEGINS (Confd) 

2.2 Presentation of Models (Confd) 

the movement of some fuel control actuator. 

Lewis, Lecourt, and Scoville used a variation of the 
ideal diesel engine in Reference (16). The torque developed 
within the engine was designated Q and the friction torque 
^fe* output torque is the difference between the two 

torques. Q^, was postulated as linearly related to fuel 
rack position, z, i.e., 

Q DE = Kz 

and the friction torque was taken as linearly related to engine 
speed , 

®fe = A 0 +A 1' N * 

It can be seen that this model does not assume constant 
mechanical efficiency. 

In Reference (17), Garvey presents a linearized different 
tial equation describing the motion of a diesel engine. He 
assumed the output torque developed by the engine to be pro- 
portional to governor position, i.e., 

T = Cz 

T - change in torque 
z - governor position 
c - governor gain 

Parker uses the same model with a delay (18), 

T = ce z 
T = Torque 
c = governor gain 
z = governor position 
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2. DIESEL ENGINE (Cont'd) 

2.2 Presentation of Models (Cont'd) 
s = differential operator 
% - tine delay 
2*3 Low Speed Diesel Model 

The ideal diesel engine model in conjunction with a delay 
in the form of an exponential lag vias chosen as the method 
to simulate a Burmeister & Wain 7K98FF low speed diesel. 

In Reference (19), Meulengracht presents the results from 
test-bed trials of a 7K98FF . From these results, and assuming 
the ideal diesel characteristics, one can get a brakehorse- 
power-speed curve and from this curve, derive the flat torque 
speed characteristics for constant fuel settings. These 
curves are all included in Appendix B. 

2.4 Medium Speed Diesel Model 

Reference (20) describes the development of high per- 
formance medium speed diesel, the Mirrless KV Major 12. 

Since the fuel map for this diesel is included in this article, 
it was selected as the diesel to model. From the fuel map, 
and using the elementary relations presented earlier, one 
can derive the engine's torque speed characteristics for con- 
stant fuel settings* In the process of drawing these charac- 
teristics, it was observed that torque could be approximated 
quite well as being proporti cnal to fuel setting. Thus 
the expression, 

T = Kx (torque) 
x = fuel setting 
K - proportionality constant 

v?a 3 used in the simulations. 



2. DI ESEL ERGINE (Cont’d) 
2.5 Delays 
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As in the gas turbine, the diesel engine does not 
instantaneously produce a change in torque with an increase 
in fuel rate. For a four-stroke cycle multi-cylinder diesel 
engine, Parker recommends a value of when N equals engine 
rpm; this is an approximation for diesel engines with two 
or more cylinders firing per revolution. (21) For the more 
general case. Reference (16) recommends 
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designed engine speed 

p : number of cylinders firing per revolution 

Discussion 



The torque-fuel setting characteristics were found to 
be linear for the medium speed -diesel but not quite linear 
for the low speed engine. The characteristics for the low- 
speed diesel were therefore stored in the computer. For each 
diesel, a change in fuel setting results in a change in torque. 
This torque was tempered by a delay in the form at a decaying 
exponential. This method is much simpler than that employed 
for the gas turbine models since diesel torque is independent 
of speed in the engine’s operating range. 

For ship propulsion studies being considered here, it was 
not considered necessary to seek a model which would produce 
the oscillatory torque that is characteristic of internal com- 
bustion engines (cf. Reference 16). 
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3. CLUTCHES AND HYDRODYNAMIC TORQUE CONVERTERS 

3.1 Introduction 

The selection of a fixed pitch propeller and a non- 
reversing prime mover will necessarily mean the incorporation 
of a clutch or torque converter in the drive train. Gas 
turbines, for example, are limited to undirectional movement; 
and some diesels, on the other hand, are designed for operation 
in two directions. It will be observed later that clutch 
engagement and disengagement times play a critical part in 
the emergency stopping of a ship. The selection of a proper 
clutch for a drive train thus is no trivial matter. In the 
same vein, the selection of a suitable clutch model to simu- 
late the emergency stopping of a ship is important in order 
to avoid erroneous results. 

3.2 Falk Karine Atrflcx Clutch 

Figure 3-1 shows the principal part of the Airflex clutch. 
The clutch is pneumatically operated drum-type dry friction 
clutch consisting of an annular rubber/f abric air gland bonded 
to an outer, flanged steel rim and lined with friction shoes 
which engage an inner iron drum when the gland is inflated. 

A newer air-cooled version utilizes extruded aluminum shoe 
blocks with axial air passages, to which the friction brake 
block liners are riveted to insulate the rubber gland from 
the heat generated at the friction surface (cf. Reference 22). 

Rubis was able to obtain from the Fall: Corporation the 
torque expressions for the clutches .(11) The torque capacity 
vras found to be directly related to the number of brake liners 
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FALK MARINE AIHFLEX CLUTCH 






Fig. ( 3 - 1 ) 
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3 . CLUTCHES AND HYDROEffMAMIC TORQUE CONVERTERS (Cont'd) 

3*2 Falk Karine Airflcx Clutch (Cont*d) 

and the net air pressure used to engage the clutch. For a 

35-inch-diameter clutch with eight liners, the maximum clutch 

torque capability, Q , is given by the expression 

c 

Q c = 2275 P net (lb-ft) 

P n( £ = net air pressure activating the clutch (psi) 

One of the important features to consider in this clutch 
is that the inward expansion of the clutch glands during 
engagement is opposed by the centrifugal force of the rotating 
clutch assembly. The retarding force has been transformed 
to an equivalent air pressure, pc, and is expressed by the 
relation 

pc = 5 + (5.9 x 10" J )Ng 2 (psi) 

Kg - speed of the frame housing the clutch glands (x'pm) 
Rubis states that typical inflation and deflation rates are 
5 psi per second and 30 psi per second respectively; and a 
typical supply pressure is 150 psi. 

3*3 SQ500 Synchroclutch (Dental Clutch Type) 

The SQ500 Synchroclutch is manufactured by the Philadelphia 
Gear Corporation. literature on the description and operation 
was furnished by the manufacturer and is included in Appendix C. 
This particular clutch is recommended for gas turbines in the 
25000 HP at 36OO rpm range (e.g. FT^A-2 & LK2500). For a 
typical application, the torque history for the synchronizing 
(friction) clutch may be approximated by: 



3 
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. CLPTCHES Aim HYDRODIMAMIC TORQUE CONVERTERS (Cont’d) 

3*3 SQ500 Synchroclutch (Dental Clutch Type) (Cont*d) 

T = T^l-©* 4 ^ 26 ) 

T = torque at time t (ft- lb) 
t = time (sec) 

T 1 = 11,600 (ft-lb) 

Both and the time constant may be changed to suit specific 
operating conditions. (23) Other information and details on 
this clutch are included in Appendix C. 

Another clutch which operates in a similar manner to the 
SQ500 is the BLH-Dynetic synchronizing clutch manufactured 
by Koppers Co, Inc. Although no torque-time history expression 
was provided by the manufacturer, information was provided 
on the operation and design of the clutch. This is included 
in Appendix C for information purposes only. 

3.4 Hydraulic Transmissions 

One of the problems inherent in a plant driven by a diesel 
or gas turbine prime mover with a fixed pitch propeller that 
is not found in a comparable steam plant is that the former 
prime movers have minimum speeds at which they can be operated. 
The Coast Guard Ship Hamilton (WPG- 715)* for example, is powered 
by two FT4-A-6 gas turbines and two Fairbanks Morse Model 
38TD-1/8 turbo-supercharged opposed-piston marine diesel engines 
driving two controllable pitch propellers. With the propeller 
at full pitch, and the diesel at idle, the vessel will attain 
a speed of approximately 12 knots, and 18 knots with full pitch 
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3 • CLUTCHES AMD OTPRODTHAMIC TORQUE COWERTBRS (Cont*d) 

3*4 Hydraulic Transmissions (Cont*d) 

and turbine at idle. (24) Richardson describes a method 
whereby the Falk clutch can be partially inflated thereby 
permitting continuous operation or maneuvering at shaft speed 
below that corresponding to minimum turbine (or diesel) idle. (25) 
For some applications, however, the heat rate may be prohibitive 
for continually doing this, and its effect on the life of the 
clutch is not clear. 

In Reference (26), Rowen describes a hydraulic trans- 
mission system that permits propeller speeds in the low rpas 
with a diesel or gas turbine. The system has a hydraulic 
reversing gear which incorporates a fluid coupling for ahead 
rotation of the propeller and a reversing torque converter for 
astern operation. In addition, a positive drive clutch is 
included for normal operation in the ahead direction at sea. 
Figure 3-2 is a cross section of the hydraulic transmission 
system. The torque characteristics on a "per unit" basis of 
the two hydraulic elements are presented in figure 3-3** 

Speed ratio, as defined by Rowen, is the ratio of the output 
speed to the input speed (positive for ahead rotation, negative 
for astern rotation). He assumed that the hydraulic element 
efficiencies were constant over the speed ratio range, and 
defined the output torque by the relation, 

T = Q - x H? m x F 
c of pm 



* "Per Unit" is defined as the actual value divided by the design point value 
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CLUTCHES AND HYDRODYNAMIC TORQUE CONVERTERS (Conf d) 

3*^ Hydraulic Transmissions (Cont*d) 

T = converter or coupling output torque 
Q q ^. = oil flow to converter or coupling 
N = input speed of prime mover 
F = a function of the speed ratio (Fig* 3-3) 

Figure 3-^ is a block diagram of the hydraulic transmission 
simulation used by Rowe niello numbers were given by him for 
the various parameters involved in the simulation. 

3*5 Discussion 

The Falk clutch and the Philadelphia Gear synchroclutch 
were selected for inclusion in the simulation because of their 
simplicity and because of their wide spread use in ship pro- 
pulsion trains. Since no expression for clutch disengagement 
was provided by Philadelphia Gear „ it was assumed that the 
clutch torque decreased as an exponential decay with a time 
constant of 5 seconds* 
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4 « PROPELLER MODELS 

4.1 Introduction 

The propeller simulation proved to be the most challenging 
aspect of this study. An excellent study of propellers under 
conditions of reversing was written by Miniovich and is listed 
as Reference (28). Miniovich questions the validity of using 
propeller characteristics measured under equilibrium conditions 
in the simulation of backing and reversing of ships. To answer 
this, he performed reversing experiments with models in a towing 
tank and measured the propeller characteristics. These char- 
acteristics Mere found to be in agreement with those measured 
in equilibrium or quasi-stationary tests. The only notable 
discrepancy appeared to lio in the zero thrust region. From 
his equations of motion, Miniovich also calculated the values 
of basic variables that characterize the reversing process 
(e.g., time required to stop the ship, time required to stop 
the propeller, and maximum values of thrust and torque) and 
found these to differ from experimental data by a maximum 
of 10 percent. It thus appears that one can use the steady 
state propeller characteristics in a simulation and get a 
reasonable accuracy. 

4.2 Propeller Characteristics 

The traditional definition of the propeller advance 
coefficient is the familiar expression, 




ND 



V. = propeller speed of advance (ft/sec) 

A 
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4. PROPELIER MODELS (Cont’d) 

4.2 Propeller Characteristics (Confd) 

N = propeller rate of rotation (rev/sec) 

D = propeller diameter (ft) 

This coefficient can be derived from a dimensional analysis 

for propellers (cf. Reference 29) • It is obvious that relating 

fixed pitch propeller coefficients to J might prove to be 

unwieldy for ship propulsion studies since J becomes unbounded 

as the shaft goes from positive rotation (providing thrust 

to move the ship ahead) to negative rotation (providing thrust 

to move the ship astern). Miniovich solved this dilemma by 

defining the propeller characteristics in terms of both J 

HD 



and its inverse, viz. 



* V A . The result is two graphs of 



propeller characteristics instead of one; when either one 
of the coefficients tends to large values he changes to the 
alternate definition. 

At the Ninth International Towing Tank Conference (1957 ) t 
Stephanson demonstrated how the J and the l/j propeller 
characteristics could be combined with a new term, called 
the modified advance coefficient defined by, 




where all the terms are defined as before. (30) (Miniovich 
ascribes a similar expression to Kalmakov in 1940). It can 
be seen that J’ will not become unbounded as N or V. approach 



zero values 
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4. PROPELIiCR MODELS ( Coni' d) 

4.2 Propeller Characteristics (Cont'd) 

In Reference (31) Baker and Patterson discuss Miniovich's 
work and present his propeller characteristics using the modi- 
fied advance coefficient. Other forms of advance coefficients 
are offered by the two authors to make the task of simulating 
a propeller easier than as performed by Miniovich. 

4.3 Modeling the Wageningen B-Screw Series 

In 1969» Van Laimeren, Van Manen, and Oosterveld presented 
a new technique of reproducing the propeller characteristics 
of the Wageningen B-Screw Series. (32) The non dimen si onalized 
thrust and torque coefficients arc presented in the form of 
a Fourier series. The series is given as a function of the 
hydrodynamic pitch angle, defined by, 




where V^jU, and D are as previously defined. The nondimen- 
sionalized thrust coefficient, C,_* , and torque coefficient, 
C q *» are then given by 



A(K) cos pK + B(K) sin 



10 C q * - M 



A(K) COS ySK + B(K) sin 



K=0 



Once the coefficients are obtained, thrust, T, in pounds, 



and torque, Q, in ft- lbs can be calculated using. 



T 
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4 . PROP EU.ER MODELS (Cont* d) 

4.3 Modeling the VJagenlngen B-Screw Series (Cont'd) 



p - fluid density 

V^, N, D as previously defined 

The factor "10” in the expression for nondimensionalized 

torque coefficient apparently through omraission, does not 

appear in two different published versions of Reference (32). 

Simulations were performed without this factor and resulted 

in excess propeller torques (10 times that predicted from 

hand calculations). Figure (40) of Reference ( 32 ) was 

reproduced by the IBM 1130 computer without the factor of 

10 in the expression. Examination of Figure (40) revealed 

that the ordinate is labeled C * and -10 C *» the minus sign 

* Q 

being introduced to avoid overlapping of curves. It was thus 
deduced that the expression for should include a factor of 10 . 

The values of the coefficients A(K) and B(K) for a series 
of 20 terras s i.e. K = 1-20, are presented in Reference (32). 

These coefficients are presented for a wide range of B-Screw 
ser5.es. Figure 4-1 illustrates how closely a series of 20 
terras matches the measured characteristics of a B 4-70 screw 
with a pitch ratio of 1.0. Figure 4-2, furthermore, illustrates 
how well a series of 10 terms approximates the measured char- 
acteristics for the same propeller. 




where 
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4. PROPELLER MODELS (Cont'd) 

rn*mnmm t * » m • ■ ^ w i nwrrm 

4.3 Modeling the Wageningen R-Screw Series (Cont'd) 



It was probably noticed that the definition of p does 
not avert the problems previously associated with the tradi- 
tional definition of J, i.e», dividing when n=o, It would 
seem more appropriate to use 



or 



f = arcsine (7^ + (0.7 ^ HD) 2 ")'*' 

0*7 HD 



- arccosine (V^ + (0*7 ^ND) y 



instead of the arctangent definition of . The IBM 

scientific subroutine library on the IK-1 1130 digital computer 

that was used did not have subroutines for calculating arccosine 

or arcsine. The arctangent subroutine had to be used. It 

should be also noted that the IBM subroutine for arctangent 

generates values only between - 'if to * 7L It is the 

2 2 . 

writer's opinion that if a numerical overflow should occur 
in the argument for the arctangent subroutine, that this sub- 
routine will yield either ^ f 2 or - ^/2 depending on the sign 
of the overflow. No numerical overflow problems wore experienced 
using the arctangent subroutine. 

4.4 C ontrollable Reversible Pitch Propellers (CRPP) 

A simple model of a CRPP was presented by Rouen in 
Reference (26). The "per unit" equation for thrust is given 
by 

2 

Thrust = (pitch x rpra) 



3 c. 



4 . PROPELLER MODELS (Cont'd) 



4.4 Controllable Reversible Pitch Propellers (CRPP) (Cont«d) 

and the "per unit" equation for torque is given by 

Torque = [o.J x (rpm) 2 ] + (p»7 + (rpm) 2 x (pitch)^J 

In Reference (33) » P.ead presents a study of gas turbines 

driving a CRPP, and includes the torque and thrust characteristics 

for the propeller. Apart from the characteristics, all that 

is known about the propeller is that it has 4 blades and has 

2 

a polar moment of inertia of 7050 lb-ft-sec . Wendel and 
Dunne apparently found a small error in the va?.ues of modified 
torque coefficient given by Read and applied a correction to 
the torque characteristic . (34) The modified advance coeffi- 
cient, J’, is as previously defined, viz., 

J* = V A 



(V^ + N 2 B 2 y 

The modified thrust coefficient is defined by 

T 

v » _ w * — ’ VT” 

*T f D (V A + H D ) 



and the modified torque coefficient is defined by 

V - 

f I^(V A 2 + n 2 d 2 ) 

where T = thrust in lbs 
Q = torque ft- lbs 

, V A , N, D are as previously defined 
(Note the absence of^ in the denominator in the expressions 
for K t ’ and K *; ft is included in the similar expressions 

1 V 
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4. PROPELLER MODELS (Confd) 

4.4 Controllable Reversible Pitch Propellers (CRPP) (Cont*d) 
given by Van Larameren . Baker and Patterson discuss the 
" r question" in Reference (31)). 

Figures (4-3) and (4-4) are the CRPP characteristics, and 
include the Uendel-Dunno correction. 

4.5 Pitch Control System Model 

The only analytical description of a CRPP pitch control 
system found was that described by Read in Reference (33) « 

Pitch, R, and angular rotation,G, are related by the expression 
R = Tan 0 

0^55 

The pitch-changing cylinder has a total capacity of 40 gallons, 
and the blade angle varies very nearly linearly with ram stroke, 
such that 

x £q (Rad/sec) 
dt 57.3 m 

Fo = oil flow in gal/ sec 

By manipulating the two equations, one can arrive at the 
expression 

dR = t 0.0528 + 0.0109 R 2 F n 
cil' 0 

There are two pumps supplying the oil flow, Fq. One is 

motor driven and the other is shaft driven. The motor 

driven pump has a maximum capacity of 60 gal/rain and the 

shaft driven pump has a maximum capacity of 100 gal/ min. 

The pumps are pressure controlled such that any alteration 

in pitch demand causes the motor-driven pump to provide 60 
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4 . PROPELLER MODELS (Cont»d) 

4.5 Pitch Control System Model (Cont*d) 

gal/ rain and the shaft-driven pump to provide its maximum 
output at the shaft speed in question* The result is that 
during a pitch alteration, 

F q = (1 + 0.435 N) (gal/ sec) 

N = shaft speed (rev/sec) 

4.6 Discussion 

The fixed pitch propeller model selected was the Fourier 
representation of propeller characteristics as presented in 
Reference (32). Only the first ten coefficients were used 
in the representation. The resistance characteristics of the 
vessel, and the propeller diameter are two important factors 
in selecting the proper propeller. Since the characteristics 
for the B-4 series with an expanded area ratio of .7 were 
readily available, it was decided to use this type of propeller. 
A plot of K^J indicated a pitch-diameter ratio of 1.0 would 
give the maximum efficiency; hence this propeller was chosen. 

Read’s controllable reversible pitch propeller (CRPP) 
was chosen. The propeller characteristics as modified in 
Reference (34) were redrawn for constant values of modified 
advanced coefficient (J*)« The independent variable was 
taken as pitch diameter (PD) ratio and the dependent variables 
were thrust and torque coefficients. Figure (4-3) shows the 
thrust coefficient vs. PD ratio for J* = -0.2. This curve 
is fairly typical of the others. The linear approximations 
are indicated. Figure (4-4) shows the torque coefficients 
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4. PROPELLER MODELS (Cont*d) 

4.6 Discussion (Cont*d) 

vs PD ratio for J* = -0.2 with the linear approximations 
used. This curve is also fairly typical of the other 
curves of torque coefficients for constant J* and varying 
PD ratio. 

Read*s pumping system as described earlier was also 
used in conjunction with the CRPP. 



5 . SHIP CHARACTERIS TICS 

5.1 Wake Fraction and Thrust Deduction Factor 

Wake fraction, w, determines the velocity the propeller 
is advancing at in relation to the ship velocity, 

V A = (1-w) V 

= propeller speed of advance 
V = ship speed 
w = wake fraction 

Thrust deduction factor, t, indicates how much propeller 
thrust is reduced due to hull interaction, 

T 1 ■ T m, <«> 

T = open water propeller thrust 
T^ = actual thrust 
t = thrust deduction factor 

(For a physical explanation of "t" and "w", and other useful 
relationships involving these two terras, see Reference (29) , 
Chapter VII, Section 13). 

The implication in the two above definitions is that 
they are steady state quantities. Thus at equilibrium, the 
actual thrust equals ship resistance. Since t«l., the 
definition shows how much more thrust is required over resis- 
tance to keep ship speed constant. 

Wake fraction and thrust deduction factor are not, as 
one would expect, constant over a range of speed. The variation 
in t and w depend on the degree of interaction between propeller 
and hull. Harvald performed several towing tank tests to 
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5. SHIP CHARACTERISTICS (Cont’d) 

5«1 Wake Fraction and Thrust Deduction Factor (Cont’d) 

measure ”t" and “w" under transient conditions for two 
different hull shapes. (35) He performed his towing tank 
tests under conditions of constant model ship velocity and 
varying shaft rpm and vice versa. Although his work pro- 
vided good insight into the behavior of "t" and "w" under these 
two conditions, it was difficult to assess how his results 
could be applied for conditions of varying shaft rpm and ship 
velocity. How did others treat this phenomena? In the solu- 
tion of his backing equations, Minicvieh assumed both factors 
to be equal to zero. ( 36 ) Rubis assumed varying values of 
"t” and "w" from two curves of bell shaped form. ( 3 ?) Read 
considered wake fraction to be negligible, but formulated 
linear relations for thrust deduction factor for positive 
and negative values of propeller thrust. ( 38 ). Lewis, Lecount, 
and Scoville assumed thrust deduction factor and wake fraction 
to be a function of vessel speed only; for negative speeds 
thrust deductions and wake fractions were neglected. (39) 

Thus, it can bo seen that there is no general agreement on 
the treatment of these two factors. The feeling appears to 
be that these two factors are so small that any errors in 
their treatment will have negligible effects on the total 
simulation. For simulations performed here, constant values 
of wake fraction (w=0.04) and thrust deduction factor (t=0.02) 



were assumed 



5. SHIP CHARACTERISTICS (Cont«d) 

5*2 Added Mass of Ship and Propeller 

As with thrust deduction and viake fraction, the treat- 
ment of added mass of ship and propeller is also varied. 

The added mass of the ship is usually taken to be between 
8-10$ of the displacement. This study used 10$. Lewis, 

Lecount, and Scoville used a different approach to account 
for added mass as will be discussed later. The added mass 
of a fixed pitch propeller is usually assumed to be about 
25$ of the inertia. (cf.,e.g., Reference (40)) Twenty-five 
percent was assumed here. There is nothing fixed about this 
figure as Burrill and Robinson found the percentage addition 
should be as high as 46$ for one particular propeller. (41) 

For his CRPPj. Read proposed in Reference (33) the expression, 
I = 325 + 1050 /R/ 

2 

I = added mass due to entrained water (Ib-ft-sec ) 

R = pitch diameter ratio 

This expression was used in the simulation with the B & W 
7K98FF low speed diesel - CRPP combination. 

5.3 Ship Resistance 

Accurate representation of ship resistance is, of course, 
vital to the proper selection of propeller (s). One of the 
basic assumptions of the popular "Propeller Law" (cf. References 
(14) & (15) ) is that resistance is proportional to the square 
of velocity for ships where the wave drag is small. For those 
cases where it is desired to simulate transients on a parti- 
cular ship, it is, of course, more accurate to use that 
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5 * SHIP CHARACTERISTICS (Cent* d) 

5*3 Ship Resistance (Cont'd) 

vessel's characteristics as determined from towing tank tests. 

Lewis, recount, and Scoville included acceleration effects 
in their expression for ship resistance, 

"h =R ss + °AD M ^ 

M = mass of ship 

V = ship speed 
s 

R " R (V ) = steady. state resistance 
ss ss s 

C. r , - "added mass" coefficient 
AD 

No values for were provided by the authors. 

For simulations performed here, resistance wa3 assumed 
to be of the form, 

R = 104.8 V 2 (lbs) 

V = ship velocity (ft/soc) 

5.4 Shaft Friction 

Read used the expression 
Q f = 2890 N 1 * 2 - 5 (ft-lb) 

N = propeller speed (assumed to be rev/ sec) 
to represent shaft friction. ( 33 ) Rubis, on the other hand, 
represented shaft friction characteristics in the form of a 
"table-look-up" in his computer simulations. These characteris- 
tics were approximated in simulations performed here by the 
expressions 

Q f = 2500 (ft-lb) 0 5/N/ £ 25 
Q f =100/N/ (ft-lb) /K/ > 25 
N = propeller shaft rpm 



40 



6. RESULTS 

6*1 Introduction 

The following maneuvers were performed with the drive 
trains as indicated: 

FT4A-2 Gas Turbine - Fixed pitch propeller 

1. Coastdown from 30 knots 

2. Accelerating from 5 knots 

3. Crashback from 30 knots 

LJ2500-A Gas Turbine - controllable reversible pitch 
propeller 

1. Accelerating from 15 knots 

2* Crashback from 30 toots 

Uirrless KV Ua.jor 12 - Fixed Pitch propeller 

1. Crashback from 9 toots 

Burmeister & Vfain 7K9&FF - Controllable reversible pitch 
propeller 

1. Crashback from 12 knots 

6.2 Coastdown from 30 toots (FT4A-2) 

At zero time, fuel rate was decreased from 9500 lb/hr to 
1600 Ib/hr exponentially with a 1-second time constant. 
Results of this simulation are presented in figures (6-1) - 
(6-4). 

6.3 Accelerating from 5 toots (FT4A-2) 

At zero time, fuel rate was increased exponentially with 
a 1-second time constant from 2500 LB/hr to 12000 LB/hr. 
Figures (6-5) - (6-11) are plots of the results. A copy of 
the computer program used is included in Appendix F. 
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6. RESULTS (Cont'd) 

6.4- Crashback from 30 knots (FT4A-2) 

This maneuver was accomplished in three phases: 

1. Coastdown* - opening clutch 

2. engaging reverse clutch 

3* stopping the ship 

Results of this simulation are plotted in figures (6-12) - 
(6-17). Time for forward clutch disengagement was 4-8 seconds. 
The reverse clutch was engaged at 11.2 seconds and the ship 
stopping time was 43*4 seconds. Head reach was 1229 ft. 

A copy of the computer program used in phase 3 is included 
in Appendix F. 

6.5 Accelerating from 15 toots (IK2500-A) 

A description of the control strategy and the controller 
used is in Appendix G® Initially the ship is at 15 knots 
with a pitch ratio of 1.0 on the CRPP. Pitch was increased 
to 1.6 in 8 seconds. At that instant the controller acted 
to increase shaft rpm to 216. Results of this simulation 
are shown in figures (6-18) - (6-20). 

6.6 Crashback from 30 knots (LI-12 500 -A) 

The problem to consider here was whether or not the power 
turbine would overspeed while the CRPP goes through a zero 
value of torque coefficient. This ship is initially at a 
speed of about 30 knots with full pitch ratio (1.6). The 



*Coastdown as used hereafter implies setting fuel rate to idle 
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6. RESULTS (Cent *d) 

6 .6 Crashback from 30 knots (LH2500~A) (Cont»d) 

first simulation that was tried used the following sequence 
of events: 

a. reduce fuel to idle 

b. decrease pitch to -1.2 

c. increase fuel to 9900 Ib/hr 

In performing this simulation, it was observed that the pitch 
changing operation took a rather lengthy tine. The end result 
would have teen a time to stop in excess of 60 seconds. The 
simulation was stopped and a new sequence of events was 
attempted, 

a. simultaneously decrease fuel to idle and decrease 
pitch 

b« when pitch equals -.8, increase fuel to 9900 
Using this sequence of events, a stopping time of 44.3 
seconds resulted. Results of this simulation are plotted 
in figures (6~2&'i - (6-24). A copy of the computer program 
used in sequence (a) is included in Appendix G. 

6.7 Crashback from 9 knots (KV Major 12) 

The same sequence of events described in section 6.4 
were used here , viz . , 

1. coastdown - opening clutch 

2. engaging reverse clutch 

3. stopping the ship 

Forward clutch disengagement occurred at 11 seconds, and 
6.2 seconds later, the reverse clutch engaged. The ship 
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6. RESULTS (Confd) 

6.7 Crashback from 9 knots (KV Major 32) (C<ont c d) 

was stopped in 95 »2 seconds with a head reach of 3.208.7 feet. 
Results of this computer simulation are plotted in figures 
(6-25) - (6—2 7 ) - Copies of the programs used in all three 
sequences are included in Appendix H. 

6.8 Crashback from 12 toots (B & H 7K98FF) 

This simulation successfully illustrated the increase 

in prime mover speed as the CRPP went through a zero torque 

coefficient region. The sequence of events employed here were: 

1. simultaneously decrease fuel to idle and pitch ratio 

2 « when pitch ratio equaled -.2, increase fuel to stop 
the ship 

The initial simulation started from a pitch ratio of 1.0 and 
decreased pitch to -1.2. At this pitch ratio, the propeller 
torques were so high that the equivalent of diesel stalling 
occurred. The change in pitch ratio was modified to go from 
+1. to -1. This proved to be satisfactory. The time of pitch 
change was 3^*5 seconds. Ship stopping time was 90.5 seconds 
with a head reach 6 29.6 ft. The results of this simulation 
are plotted in figures (6-28) » (6-31). A copy of the pro- 
gram used in sequence (2) is included in Appendix I. 
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?. CONCLUSIONS 

Simple models wore found to exist that enable one to model 
ship propulsion drive trains. Some of the models had to bo refined 
by the addition of time lags, delays, etc. The combination of 
drive train components must be compatible with the characteristics 
of each drive train member. The simulations illustrated the 
importance of proper sequencing of propulsion plant control para- 
meters. Conditions of over speeding, excess torques, and high 
energy dissipation rates in clutches can be easily simulated. 

These simulations can prove to b3 immensely useful to the control 
engineer in helping him synthesize a controller that will keep a 
system within its region of safe operation. 
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8. RECOMMENDATIONS 

1* A study should be taken to determine what is the optimal 

fuel flow getting) - pitch ratio combination for a particular 
prime mover and the CBPP presented here. 

2. A model of reversing the B & W 7K93FF should be developed 
and a crashback simulation of the 7K98FF with a fixed pitch 
propeller should be performed. 

3. A model controller to keep shaft rpm nearly constant while 
changing pitch ratio on the CRPP could be synthesized for 
each one of the prime movers. 

4. A torque converter model should be developed to be used with 
the prime movers developed here, and maneuvers should be 



simulated 



REFERENCES 



Kiviat, P.J., Digital Computer Simulation: Modeling Concepts , 

RM- 5378 -PR (Rand Memorandum )~, August, 1967 ; AD 656429. 

Kurzak, K. H», et. al., Propulsion Machinery of the "Koeln" 

Class Escort Frigates With Special Consideration of Gas Turbine 
Propulsion , ASHE Publication 65-G TP-11 . 

Cottingtcn , R , V . , The Da sign o f a Digital Three-Term Controller 
as a Turbojet Engine 5 ;p»ed Governor' usin'A Digital Simulation 
Methods , National Gas Turbine Establishment (Great Britain),' 
Reports and Memoranda No. 3615 , November, 1967. 

G. S. Mueller, "Gas Turbine Simulation Using One-Diraensional 
Flow Relationships," XSE C onference on Industrial Applications 
of Dynamic Modeling , ±6^1& September , 1969 > University of 
Durham; p. 123 • 

Marscher, William F«» Lin ear Analysis of Aircraft Gas Turbine 
Control System Dynamics^ Master of Science Thesis, Department 
of Mechanical Engineering, Massachusetts Institute of Technology, 
I960. 

Arthur D. Little , Inc., First Interim Technical Report, Computer 
Slmiilati on of the Power Conversion System for Ijike-X Power plant s 
DA-49-129-Eng-542 

Saravanamuttoo , K.I.H., Fawke, A. J., "Dynamic Modelling of Gas 
Turbine Performance," IGE Conference on Indus tri al Applications 
of D ynamic Modelling , University of Durham, l5^l8 September, 1969 
P * T33 • 

Saravanamuttoo , Fawke, A.J., Simulation of Gas Turbine 

Dynamic Performance , ASMS Publication 7O-GT-23. 

Bichayev, B. P., Analog and Digital Models of Marine Gas-Turbine 
Units , Shipbuilding Press, Leningrad, 1969 , in Russian. 

Lnrye , I. Yu, "Mathematical Modeling of a Gas Turbine with the 
Use* of a Real Fuel Apparatus," Sudostroyeniye, No. 12, 1964, 
pp 24-2.6, AD 645353. 

Rubis, C. J., "Braking and Reversing Ship Dynamics," Proceedings 
Second Ship Control Systems Symposium , November 4-6, I 969 T” 
sponsored by Naval Research and Development Laboratory (Annapolis 
Mid.), p. IV-B-2. 

Bodnaruk, A., Rubis, C. J., Acceleration Performance of a Gas 
Turbine Destroyer Escort, ASMS Publication ?0-GT-77» 



86 



REFERENCES (Cont'd) 



( 13 ) 

< 1*0 

(15) 

( 16 ) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 



Rubis, C. J., Transient Re sponse of a 25,000 HP Marine Gas-Turbine 
Engine t ASMS Publication 71-GT-61. 

Woodward, J.B. ,111, Tine Diesel Engine: To Drive A Ship , Department 

of Naval Architecture and Marine Engineering Publication No, 105> 
January, 1971* University of Michigan. 

Pounder, C. C., Marine D iesel E ngines , Newnes-Butterworth, London, 
pps. 186-258. 



Lewis, J. W«, et.al., Digital Computer Simulation of Diesel-Electric 
Propulsion Systems , Presented during the Sep tembe v~2&^2 ?7 ' & 2<T~ 
Joint meeting of the Chesapeake & Hampton Roads Section of the 
Society of Naval Architects and Marine Engineers, 

Garvey, D., Feedback Control Systems Analysis , Paper No. PKCC 
71-11, Woodward Governor Company (PKCC is believed to stand for 
Prime Mover Control Conference, a meeting held annually at the 
Woodward Governor Company). 



Parker, George, G eneral Discussion and Eouatlons of 
Paper No. PKCC ?6-4 , Woodward Governor Company ( same- 
in previous reference). 



a Prime Move 
comments as 



Kenlengracht , Per V. .Hi gh-Powered Marine Diesel Engines , Paper 
presented in Vancouver^ British Columbia, March 2," 1970 to a 
combined meeting of The Institute of Marine Engineers and the 
Society of Automotive Engineers. 



Pope, J.A., Lowe , W. , "The Development of ?. Highly-Rated Medium- 
Speed Diesel Engine of 7000-9000 Horsepower for Marine Propulsion 
Transactions , Institute of Karine Engineers, Vol. 78, No. 8, 1966. 

Parker, op. clt. , p.2. 

Richardson, W. S. ? The Friction Clutch Reverse-Reduction Gears 
for the GTS, ASMS Publication 69-01-5. 

Private correspondence dated April 10, 1972, from J. K. Liu, 
General Manager, Synchrotorque Division, Philadelphia Gear 
Corp oration. 

Wennberg, P. K., "The Design of the Main Propulsion Machinery 
Plant Installed in the USCGC Hamilton (WPG-715) * " Transactions 
SNAKE, 1966. 



87 



REFERENCES (Cont'd) 



(25) Richardson, W. S. op.cit . , p. 5» 

(26) Rowen , W# I., A Comparison of S imulated Gas Turbine Ship Handling 
Characteristics with Several Different Transmission Syst ems, ASMS 
Publication 71-GT-65. 

(27) Ibid ., p. 7. 

( 28 ) Miniovlch, I. Ya. } Investigation of Hydrodynamic Characteristics 
of Screw Propeller Under Conaiticns oi Re v ersir g~ana Calculatio n 
Met hods l’or Tacking of Snips , Bureau of Ships Translation No. 597 
(translation by Roger & Roger Inc.), originally appearing in 
Transactions of the Red Banner Central Scientific Research Institute, 
Issue No. 122, published by the State Union Publishing House of 

the Shipbuilding Industry, 1958* 

(29) Comstock, John P. t ed*, P rinciples of Naval Architecture , The 
Society of Naval Architects and Marine Engineers, 1967, p. 385 * 

(30) Stephanson, T., Propeller Constraints , Ninth International Towing 
Tank Conference, 1957* 

(31) Baker, Dennis W., Patterson, Clair L., Jr., "Some Recent Develop- 
ments in Representing Propeller Characteristics," Proceedings 
Second Ship Control systems Symposium, November 4-5-6, 

pps « AP-A-1 - AP-A-34. 

(32) Van Lansnsren, W.P.A., Van Kanen, J.D., Oosterveld, M.W.C., The 
Wageningon B-Screw Series , presented at the Annual Meeting, New 
York, N.Y., November 12-24, 1969 , of the Society of Naval Architects 
and Karine Engineers. 

(33) Read, C.M.B. , "The Control of Gas Turbines Driving Controllable- 
Pitch Propellers," P roceedings of the First Ship Control Systems 
Symposium, Vo 1. I, November 15-17 > 1966, held at Annapolis, Maryland. 

(34) Wendel, Alan H., Dunne, James, F., " Dynamic Analysis and Simulation 
of Ship and Propulsion Plant Maneuvering Performance," Pro ceedings 
of the second Ship Control Systems Symposium, Vol. I, HovemberT^b, 
1969 , held at Annapolis, Maryland. 

(35) Harvald, SV. AA., 'Wake and Thrust Deduction at Extreme Propeller 
Loadings," Publications of the Swedish State Shipbuilding Experi- 
mental T ank, Not 6lf~ 1967 ™ 



( 36 ) Miniovich, I. Ya., op* cit , p. 33* 



REFERENCES (Cont'd) 



(37) Rubis, C. Jo, op, cit «, p, IV-B-4. 

(33) Read, C. Mo B., op. clt ,, p. VII-B-8. 

(3S) Lewis, et. al., op, cit, , p, 10, 

(40) - Norrbin, Nils H«» "Theory and Observations on the Use of a 

Mathematical Model for Ship Maneuvering in Deep and Confined 
Waters,." Publ ications of the Swedish Stato Shipbuilding 
Experir:entaX~Tan k , TjoTTbT** 

(41) Bur rill, L. C. , Robson, W«, "Virtual Mass and Moment of Inertia 
of Propellers, " Tra ns actions North East Coast Institution of 
Engineers and Shipbuilders, Vol. ?8, 196 1-62, pps. 325-350, 



89 



APPENDIX A 
Gas Turbine Models 

1. Lurye Model 

The following model of a marine gas turbine was proposed by 
Lurye in Reference (10) in the form of linearized differential 
equations: 

Turbocompres sor 
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APPENDIX A (Cont'd) 



1 . 



(Cont'd) 
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.... load coefficient by propeller 
cJ/^ pitch 



c)<^> 






.... amplification coefficient by 
jyy rate of compressor rotation 





turbocompressor rate of rotation 
free turbine rate of rotation 



G fj , = fuel consumption 
p - propeller pitch 
£ = time lag in combustion chamber 





M 2 

\ 



— f 

Tj_ (cv^Gy) torque moment of turbocompressor 

torque moment of free turbine 

- turbocompressor rotational inertia 
= free turbine rotational inertia 

= f j ft ) resistance moment of propeller 

= fj, O resistance moment of compressor 
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APPENDIX A (cont’d) 



2. An Extension of the Saravanamuttoo Model 

Using the technique proposed by Saravanamuttoo in References 
(7) & (8) $ the model can be extended to the case where 2-spool gas 
generator-free turbine combination: 

HP Rotor equation 






Kr-r C ,„„ (7i - Z) - 



IP Rotor equation 







Free Turbine equation 







c/X 



LP-HP Intercompressor Volume 




rt 2 d % 



v 2 c/*r 



Combustor volume 



W - W - VI „ 
HPT HPC f6 
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Interturbine Volume (IP -HP) 




% sr 



Interturbine Volume (HP-FT) 
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FT 



c/d 
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(Cont'd) 

Combustor 

T 4 ~ < ~v c < ^3 

Nomenclature 

W = mass flow rate 
T = temperature 

Cp= specific heat at constant pressure 
?L r - turbine over-all efficiency 
^ - compressor over-all efficiency 
K - gas constant 
HV= fuel heating value 
V = volume 
U = rotational speed 
P = pressure 
t = time 

I = polar moment of inertia 
Station Numbering 

1 = IP Compressor inlet 

2 = HP Compressor inlet 

3 = combustor inlot 

4 = HP turbine inlet (combustor outlet) 

5 = IP turbine inlet 

6 = FT inlet 



7 = FT outlet 
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(Cont'd) 

Subscripts 

LPC - low pressure compressor 
ILPC = high pressure compressor 
FB ~ Combustor 
KPT ~ high pressure turbine 
LPT = low pressure turbine 
FT = free turbine 
HP = HPC-HPT rotor 
LP = LPC-LPT rotor 
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APPENDIX A (Cont’d) 

3* FT4A-2 & LM2500-A Marine Ga3 Turbines 

Figure (A-l) is the pcwer-speed curve for constant SFC values. 
The torque-speed characteristics shown in figure (A~2) can be 
derived from this curve. Figure (A-3) and (A— 4 ) are the corres- 
ponding curves for the LM2500-A. Figures (A-2) and (A-4) show 
examples of the linear approximations for constant fuel rate 
values used to approximate the gas turbines* characteristics. For 
intermediate values of fuel rate, linear interpolation vras used. 

The subroutines used for the FT4A-2 and the LM2500-A are included 
in the rear of this appendix. 

A sketch of the FT4A-2 with a JFC 25-23 fuel controller is 
shown in Figure (A-5)c A controller similar to the JFC 25-28 is 
the JFC 25-7 shown in figure (A -6). This controller is typical 
of the non-linear systems engineers are often asked to model. 
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SEA LEVEL 



T AMBIENT = 59° F 

BASED ON LIQUID FUEL LHV = 18, 500 BTU/LB 
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SHAFT HORSEPOWER vs. POWER TURBINE SPEED FOR LH2500-A 






Fig. (A-3) 
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HAMILTON STANDARD 

JFC 25-7 FUEL CONTROL 
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PACE 1 

SUBROUT I ME LM2 50 ( RPM , WF , TORQ , A08 » A 1 6 , A2 , A3 » A 4 , A5 , A6 , A7 , A99 ) 
THIS SUBROUTINE CALCULATES FREE TURBINF TORQUE AS A FUNCTION OF 
RPM AND FUEL FLOW 
MINIMUM FUEL RATE IS BOO LP/HR 
MAXIMUM FUEL RATE ALLOWABLE IS 9900 LB/HR 

ENTERING ARGUMENTS FOR THIS SUBROUTINE ARE FUEL RATE (WF) 

AND POWER TURBINE RPM (RPM) OUTPUT IS POWER TURBINE TORQUE IN 
FT-LBS 

SEARCH TO BRACKET WF AND THEN PROCEED TO FUf : L SUBROUTINES 
I F ( WF .FQ. 800.) GOTO 15 
I F ( W F .LE. 1600.) GO TO 1 6 
IF(WF .LE. 2000.) G0T017 
I F ( W F .LF. 3000 . ) GOTO I 8 
I F ( W F .LE. AOOO.lGCrOLO 
I F ( W F .LE. 5000.) GOTO II 
I F ( W F .LE.6000. 1G0T012 
I F ( W F .LE. 7000 . ) GO TO i 3 
I F ( W F .LE. 9900 . ) GOTO 1 4 

15 T0RQ = -1.965*RPi v + 5500. ■ 

RETURN 

16 CALL LWF08 (RPM, A08) 

CALL LWF16 (RPM, A 16) 

T0RQ=A08+( WF-800. ) * . 00 1 2 5 » ( A 1 6- AO 8 ) 

RETURN 

17 CALL L W F 2 (RPM,A2) 

CALL LWF16 (RPM, A 16) 

T0RQ=A16+ ( WF- 1600. ) * .0025* ( A2-A16 ) 

RETURN 

18 CALL LKF3 (RPM, A3) 

CALL LWF2 (RPM,A2) 

T0RQ=A2+t WF-2000. ) * . 001* ( A3-A2 ) 

RETURN 

10 CALL LWF4 ( R P M , A A ) 

CALL L W F 3 (RPM, A3) 

T0RQ=A3+( WF-3000 . ) « .001* ( A4-A3 ) 

RETURN 

11 CALL LWF5 (RPM,A5) 

CALL LWF4 ( R P M , A A ) 

T0RQ=A4+( WF-AOOO. ) * . 00 1* ( A 5- A A ) 

RETURN 

12 CALL LWF6 (RPM.A6) 

CALL LWF5 (KPM.A5) 

T0RQ=A5+( WF-5000. ) *.001* ( A6-A5) 

RETURN 

13 CALL LWF7 (RPM,A7) 

CALL LWF6 (RPM,A6> 

T0RQ=A6+( WF-6000. ) * . 00 1 * ( A7- A6 ) 

RETURN 

1A CALL LWF99 (RPM.A99) 

CALL LWF 7 (RPM.A7) 

T0RQ=A7+ ( WF-7000. ) * . 0003AA8 * ( A99- A7 ) 

RETURN 

END 

SUBROUTINE LWF08 (RPM,A0'8) 

A08=-1.965*RPM+5500. 

RETURN 



05/03/ 
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PAGE 2 05/03/ 

END 

SUBROUTINE LWFI6 (RPM,A16) 

A16=-2.94*RPM+12500. 

RETURN 

END 

SUBROUTINE LWE2 (RPM,A2) 

A2 = -3.48*RPFi +16150. 

RETURN 

END 

SUBROUTINE LWF3 (RPM,A3) 

A 3 = - A .77*RPM +23250. 

RETURN 

END 

SUBROUTINE LWF4 (RPM,A4) 

A4=-5. 933* R PM +30800. 

RETURN 

END 

SUBROUTINE LWF5 (RPM,A5) 

A 5 = -7 . 133*RPM+39500. 

RETURN 

END 

SUBROUTINE L W F 6 (RPM.A6) 

A6=-7.5*RPM +47500. 

RETURN 

END 

SUBROUTINE LWF7 (RPM,A7) 

A7=-8.*RPM +55300. 

RETURN 

END 

SUBROUTINE LWF99 (RPM,A99) 

A99=- 10. 1 5*RPM+73600 . 

RETURN 

END 
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APPENDIX B 
Diesel Particulars 



1 . 



2 , 



Burmeister and Wain 7K9&FF 

Test-bed results of the B & W 7K9SFF are shown in Figure 
(B-l). From this figure, it was possible to derive the ideal 
diesel power-speed curves for constant fuel settings; this is 
shown in figure (B-2). The torque-speed characteristics, figure 
(B-3). are in turn derived from figure (B-2). The "FF" at the 
end of the designation of this engine denotes that this diesel is 
intended for direct drive, i.e., no reduction gears. Other 
technical data is presented below; 

continuous service rating 24500 metric BHP @ 100 rpm 
maximum continuous service rating 26600 metric BHP © 103 rpm 
Length 45 • 2" 

Breadth 16* 2" 

Height 43' 2” 

The computer program used to model the diesel is included 
at the end of this appendix, 

Mirrless KV Major 12 

A cutaway drawing and technical specifications of this engine 
is shown in figure (B-4). The "12" at the end of the designation 
indicates the number of cylinders. Figure (B-5) shows the fuel 
map for this engine. Using this fuel map and the elementary 
diesel equation discussed in section 2, it was possible to obtain 
the approximately Linear torque-fuel setting characteristics shown 
in figure (B-6). The relationship obtained from figure (B-6) is 
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APPENDIX B (Cont'd) 



2. (Cont’d) 

T = 55 0z 

z = fuel setting in percent 
T = torque in ft-lbs 



RPM 




Fit;. v£~l) 






no 




B.m.z.p fb./sq.it 
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Fig. (b~5) 



SUBROUTINE BURN E (WF,7CRQ) 

THIS SUBROUTINE COMPUTES 7K98FF DIESEL TORQUE AS A FUNCTION 

FUEL SETTING. 
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APPENDIX C 



Clutches 

Figures (C~l) and (C-2), which were provided by the Philadelphia 
Gear Company, describe the operation of quill shaft version of the firm’ 
synchroclutch. Additional useful information provided by Philadelphia 
Gear Is: 

Inertia . . Input -17^ lb-ft^ 

output -no ib-rr 

torsional stiffness*.......... Approx. 270 million in-lb/Rad 

Figures (C~3) and (C -4), which were provided by the Koppers Co., 
describe the operation of the BUI-Dynetic clutches. 



* For information purposes only 
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A FRICTION CLUTCH PACK 
B ROTARY UNION 

C. OENTAL CLUTCH ASSEMBLY 

O. ENGAGING CYLINOER FOR DENTAL 
CLUTCH 

E. THRUST BRG. 

F. SHIFT COLLAR 

G. PINS 

H. SLEEVE 

I INTERMEDIATE CEAR 
) DRIVE SPLINE 
K. OUTPUT FLANGE 
l SELF LOCK PISTON ASSY. 

M DISENCAGE CYLINDER-DENTAL 
CLUTCH 

H SPACE FOR OPTION. QUILL SHAPI 
BRAKE 

D. FRICTION APPLY CYLINDER 

P. HYDRAULIC BALANCE 4 FRICTION 
RETURN PISTON 

Q. JOURNAL 4 THRUST BEARING 
ASSEMBLY 



FIG. 2 CROSS SECTION OF TYPICAL SYNCHROCLUTCH {QUILL SHAFT VERSION) 



DESCRIPTION i 

B The new Philadelphia SYNCHROCLUTCH consists of a 
friction clutch in parallel with a positive drive dental clutch. 
The friction clutch synchronizes the driving and driven 
shafts, after which the dental clutch is engaged to provide 
a positive drive. The dental clutch features fully crowned, 
nitrided gears for maximum capacity for a given diameter. 
The friction clutch is of the multi-disc type with positive 
oil circulation for long life and heat removal. The friction 
clutch has ample capacity to absorb thermal loads imposed 
by the synchronizing cycle— even very high loads encoun- 
tered when reversing large marine gears in the 40,000 
plus HP range. In addition, dead shaft pickup is possible 
in large marine propulsion systems. This permits the prime 
mover to be run at speed in port for auxiliaries; such as 
generators, while retaining the ability to get underway 
instantly by clutching in the propeller. 

The dental clutch features a unique helical gear element 
which automatically aligns the teeth of the dental clutch 
prior to engagement. This feature is illustrated in Fig. 3. 
The clutch may bo actuated hydraulically or pneumatically, 
using low pressure air (100 psi) or oil (150 psi). The 
synchronizing action is controlled by a transistorized dif- 
ferential speed control using non-contacting magnotic 
pick-ups for speed indication and position readout. Once 
engaged, the dental clutch does not require power to 
remain in mesh, thus eliminating dependence on electrical 
or hydraulic power for continued torque transmission. 




External jacking is provided for emergency manual engage- 
ment and disengagement of the dental clutch. 

Referring to Fig. 2, the friction clutch (a) is first engaged 
to force the driving and driven shafts into synchronism. 
The clutch Is camped by admitting oil (or air) to the piston 
through the rotary union (b). After the driving and driven 
shafts have reached a speed differential of approximately 
1/2 % . the transistorized speed control signals for dental 
clutch (c), engagement by admitting oil (or air) pressure 
to the engaging side of the dental clutch piston (d). The 
piston transmits the thrust through a heavy duty thrust 
bearing (e), to the rotating shift collar (f), which is con- 
nected through large pins (g), to the sleeve (h), of the 
dentai clutch. The sleeve carries with it the intermediate 
gear (i), and meshes with the drive spline (j), on the output 
flange (k), in the manner illustrated in Fig. 2. A self-locking 
piston (I) on the dental clutch prevents piston back-off 
under torque reversal. This action is shown in Fig. 3. To 
disengage, pressure is admitted to the reverse side (m) of 
the dental clutch piston, which pulls the dental elements 
out of engagement as shown in Fig. 3. The friction clutch 
may then be de energized to permit the driving and driven 
shafts to rotate relative to each other. Where desired, a 
quill shaft brake (n) may be incorporated into the clutch 
to stop the quill shaft from rotating, due to residual drag 
torque within the friction clutch. This feature requires no 
additional length. 



Fig* (C-l) 
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1. DENTAL CLUTCH DISENGAGED 




3. SLIDING SLEEVE HELICAL GEAfl ROTATES 
INTERMEDIATE GEAR INTO ALIGNMENT WITH 
DRIVE SPLINE TEETH 




2. INTERMEDIATE GEAR & DRIVE SPLINE GEAR 
BUTTING 




~~ 1 



J 




1. 


SPEED DIFFERENTIAL (INPUT/OUTPUT). 


UP TO 100% OF RATED SPEED. 


2. 


ACTUATION FLUID. 


OIL @ 150 PSI. (SPECIAL-AIR @ 100 PSI.) 


3. 


FRICTION MATERIAL. 


PGC. FS001 GRADE RC1. 


4. 


STEEL DISC MATERIAL. 


HARDENED STEEL. HIGH CARBON. GROUND. 


5. 


ROTATING PARTS. 


FORGED STEEL, HEAT TREATED AND 
DYNAMICALLY BALANCED. 


6. 


DENTAL CLUTCH MEMBERS. 


NITRIDED STEEL. 


7. 


DENTAL CLUTCH TEETH 


CROWNED, CHAMFERED ON ENTRY SIDE. 


8. 


ROTARY UNION. 


EXTERNAL— UNIT CONSTRUCTION. 


9 


ACCESSIBILITY. 


ALL CLUTCH MEMBERS REMOVABLE WITH- 
OUT REMOVING HUB OR QUILL SHAFT. 


10 


DIRECTION OF ROTATION. 


CW OR CCW. 


ii. 


PERMISSIBLE AXIAL MOTION OF OUTPUT 
FLANGE CONNECTION. 


±1/8" STANDARD. 


12. 


PERMISSIBLE RADIAL RUNOUT AT OUTPUT 
FLANGE. 


.015 TIR ON 500 SIZE AND LARGER. 
.010 TIR ON 120 SIZE. 


13. 


ELECTRICAL POWER REQUIREMENTS FOR 
CONTROLS. 


115V, AC-50/60 HZ, 50 WATTS. 


14. 


CLUTCH ENGAGEABLE AT STANDSTILL. 


YES. 


15. 


EMERGENCY MANUAL CLUTCH OPERATION. 


YES. 



OPBKATION AND 



SPECIFICATION OF SYNCJ1ROCLUTC11 



Fig. (C-2) 
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DISENGAGED POSITION 




FUNDAMENTAL DESIGN 

All BLH-Dynetic Clutches are of the same basic 
design, modified as necessary to meet a wide vari- 
ety of application conditions. 

The forces developed during engagement are 
completely self-contained, and produce no thrust 
loads on the shafts or clutch housing. After engage- 
ment, power is transmitted through hardened and 
nitrided coupling teeth in positive mesh, with no 
loss due to slippage. Hydraulic pressure is required 
only during the actual engaging or disengaging 
cycle. The parts are held in their final positions by 
spring loaded detents, or mechanical locks, as 
necessary. 



MAJOR COMPONENTS 

[ Operating Cylinder— to provide the axial force 
necessary to engage or disengage the clutch. Actu- 
ation is usually hydraulic using lubricating oil as the 
fluid, but pneumatic operation is also available. 

Friction Discs— to develop friction torque neces- 
sary to synchronize the shaft speeds. These discs 
are cooled by a continuous flow of lubricating oil. 

Transition Torque Control— spring loaded pins 
key with lugs in the load end hub to transmit torque 
during the transition from friction to gear tooth 
drive, and to align gear teeth prior to engagement. 

^ Gear Coupling— to transmit power after engage- 
ment is complete and accommodate shaft misalign- 
ment, if present. 



DESCRIPTION OF BLH DYNETIC CLUTCH 



Fig. (C-3) 
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The heart of the synchronizing mechanism is the 
cluster of synchronizing pins (SP) which cooperate 
with the holes in the power end hub (PEH) through 
which they project to control engagement as indi- 
cated by the diagrams at the left. 

As engagement begins the actuating piston moves 
the power end hub and its connected parts toward 
the load end hub. Through detent action, an initial 
friction torque is developed which swings the small 
diameters of the synchronizing pins against the 
sides of their respective holes in the power end hub, 
as in Figure A. 

Continued pressure brings the cam surfaces of 
the synchronizing pins into contact with those of 
the power end hub, thereby transmitting full piston 
pressure to the friction disc pack as in Figure B. 
The friction torque so developed acts to synchronize 
shaft speeds. The geometry of the assembly is such 
that the parts remain in the position shown as long 
as slipping occurs in the friction discs. 

Slipping stops at the friction discs when the shaft 
speeds synchronize. This unbalances the forces and 
permits the power end hub to advance over the syn- 
chronizing pins as in Figure C. At the same time, 
spring loaded pins are moved into spaces between 
lugs in the end of the load end hub, acting as loose 
keys to control torque until the coupling teeth 
engage. ■ 

As the synchronizing pins become aligned with 
their mating holes as in Figure D, the cam surfaces 
separate, removing the axial force on the friction 
disc pack. This permits the spring loaded pins to 
contact the sides of their lugs, maintaining syn- 
chronism of shaft speeds and aligning the coupling 
teeth for engagement as the piston completes its 
travel. 

With completion of engagement the hydraulic 
pressure is cut off, and the assembly transmits 
power as a simple gear-type coupling. 

During disengagement the piston moves the parts 
a’way from the load end hub, disengaging the spring- 
loaded pins from their lugs before the coupling 
teeth disengage, and securing the assembly in the 
final position to prevent dynamic unbalance. 



OPERATION OF IJLII. DYNSI'IC CLUTCH 



Fig. (C-4) 
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APPENDIX D 

Propeller Subroutines 

This appendix contains tho computer programs used to model the 
B4-?0 (P/d = X.) propeller and Read's CRPP, both of which are discussed 
in Section 4. Included in this appendix is the program described in 
Section that computes shaft friction. 
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**C0S(8. *8)-.0096650*SIN(8.*B)-. 03 3291* CCS ( 9 .* B ) + . 0431 9*S IM9.*B) + 
50 12 3.11* CCS ( 10.*B)+.012453*S INI 10. *6) 

CC=.1*CC 

RETURN 
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SUBROUTINE CRPP ( J , C 0 , C T , AP I , AP 8 , BP l , BP8 , AMO , AM l , BM8 , BM 1 , AMA , BMA , A 
2M2 t BM2 , AMO , BMO , APA , BPA , PD ) 

C CRPP COMPUTES MODIFIED THRUST AND TORQUF COEFFICIENTS FOR A. GIVEN 

C VALUF OF MODIFIED ADVANCE COEFFICIENT AND PITCH DIAME1ER RATIO. 

C 

C PROPFLLOR DATA IS TAKEN FROM DYNAMIC ANALYSIS AND SIMULATION 

C OF SHIP AND PROPULSION PLANT MANUVER I NG PERFORMANCE BY WENOEL 

C AND DUNNE APPEARING IN PROCEEDINGS OF THE SECOND SHIP CONTROL 

C SYSTEMS SYMPOSIUM. 

C 

c 

C MODIFIED ADVANCE COEFF J=VA/SQRT(VA**2+(N*D)**2) 

C MODIFIED THRUST COEFF . . . CT=T/ ( RHO« ( D**2 J * ( VA« *2+ ( N*D ) **2 ) ) 

C MODIFIED TORQUE COEFF. . .CQ=Q/ (RHO* ( D**3 )* (VA**2+ ( N*D) **2 ) > 

C 

C 

C RHO=DENS I TY OF WATER 

C D= PROPEL LOR DIAMFTER 

C V A = PROP E L LFJR SPEED OF ADVANCE 

C N= PROPF L LOR ROTATION RATE 

C T = T H R U S T 

C Q=TORQUE 

C 
C 

C PITCH DIAMTER RATIOS ARE LIMITED BETWEEN +1.6 TO -1.2 

C ENTERING ARGUMENTS ARE MODIFIED ADVANCE COEFFICIENT (J) AND 

C PITCH DIAMTEFR RATIO (PD) 

C 

C OUTPUTS ARE CQ AND CT 

REAL J 

IF(J.LE.-.8)G0r010 
I F ( J.LE.-.A )G0T01 l 
IF(J.LE.-.2)G0T012 
I F ( J.LE.O. JG0T013 
I F ( J.LE. .A JGOTOIA 
IFtJ.LE. . 8 ) GOTO 1 5 
CALL PLUS l (PD.AP1»BP1) 

CALL PLUS3 (PD, APB, BPS) 

F= ( J-.8 ) *5. 

CQ=AP8+F*( AP1-AP8) 

CT=BP8+F» ( BP1-BP8 ) 

RETURN 

10 CALL MIONE (PD,AMi,BMl) 

CALL KIN8 ( P D , A M 8 , B M 8 ) 

F= ( J+l . ) *5. 

C0 = AM1 + F* ( A M 8 — A M 1 ) 

CT=BM1+F* ( BM8-0M1 ) 

RETURN 

11 CALL MINA ( PD , AMA , BMA ) 

CALL MIN8 (P0,AMR f BM8) 

F = ( J + . 8 ) »5. 

C0 = AM8 + F* ( A Mi A- AM B ) 

CT=BMR+F* ( BMA-BM8 ) 

RETURN 

12 CALL M1N2 (PD,AM2,RM2) 

CALL MINA ( PD , AM A , BMA ) 
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F=( J+.4 )»5. 

CQ=AM44 F* ( A i' ' 2 - A f 4 ) 

C T = BM4 + F* ( BM2-BM4 ) 

RETURN 

13 CALL ZERO ( P D » A M 0 » B M 0 ) 

CALL MIN2 (PD,Ay.2,BM2) 

F=( J+.2 ) *5. 

CQ=AK2+F* ( AMO-AM2 ) 

CT = BM2 + F*( 800-802 ) 

RETURN 

1 A CALL PLUS4 <PD,AP4,BP4) 

CALL ZERO (PD,AMO,BMO) 

F = J * 2 . 5 

CO=AMO+F* ( AP4-AMO ) 

CT=BMO+F*( 8P4 -BOO ) 

RETURN 

15 CALL PLUSS (PD,AP8,BP8) 

CALL PLUS4 (PD,AP4,BP4) 

F=( J-.4)*2.5 
C0=AP4 + F* ( AP8-AP4 ) 

CT = 3P4 + F* ( BP8-BP4 ) 

RETURN 

END 

SUBROUTINE PLUS1 (PD,APl,BPi) 
IFtPD.GE. . 63 ) GOTO 10 
AP1=-.477*PD-.015 
G0T011 

10 AP1= .00516*PD-.049 

11 BP 1 = - . 4 
RETURN 
END 

SUBROUTINE PLUS8 (PD,AP8,BP8) 
I F ( P D . G E . *35) GOTO 1 0 
AP 0 = - . 06* PD 
GOTO 1 1 

10 AP8=.026*PD-.Q3 

11 IF(PD.LT.-.4) GOTO 1 2 
IFtPD.GT. .4) GOT 0 1 3 
BP8=.0R5*PD-.275 
RETURN 

12 B P 8 = . 165*PD-.243 
RETURN 

13 BP8= . 24 3* PD- .3382 
RETURN 

END 

SUBROUTINE PLUS4 (PD,AP4,BP4) 
IFtPD.GT. * 4 ) GOTO 10 
IFIPD.LT. - . 4 ) GOT Oil 
A P4 = . 0 1 5 
GOTO 12 

10 AP4=.0848*PD-.019 
GOTO 12 

11 AP4=-. 123*PU- .0325 

12 IFtPD.LT. - . 4 ) GOT 013 

I F ( PD .GT. .4) G0T014 
BP4 = . 24 5 * PD~ . 9 
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RETURN 

13 BP4=.37»PD-.04 
RETURN 

14 B P4= . 4 l 5 * P D- . 150 
RETURN 

END 

SUBROUTINE ZERO (PD.AMOtBMO) 

IFIPD.LT. ~.4)GOT010 
IFIPD.GT. .4) GOTO 11 
AMO =.0125 
GOTO 12 

10 AM0=-. 125*PD-.0375 
GOTO 1 2 

11 AMO = . 1092*PD- .03 1 1 

12 IFIPD.LT. -.625IG0T013 
IFIPD.GT. .425) GOTO 1 4 
0 M 0 = . 3 1 0 * P D + . 0 2 0 
RETURN 

13 8MO=.53*PO+. 1575 
RETURN 

14 BM0=.53*PD~.O732 
RETURN 

END 

SUBROUTINE MIN2 I PD , 4.M2 , BM2 ) 

IFIPD .GT. .045) GOTO 10 
IFIPD .IT. - . 04 7 5 ) GOTO 1 1 
AM2 = .0125 
GOTO 1 2 

10 AM2=.123*l>D-.043 
G0TO12 

11 AM 2=-. l32*PD-.05 
G0T012 

12 IFIPD .LT. - . 625 ) GOTO 1 3 
IFIPD .GT. . 425 ) GOTO 14 
BM2=.310*PD+.02 
RETURN 

13 BM2=.53^PD+. 1575 
RETURN 

14 BM2= . 53*P[J- .0732 
RETURN 

END 

SUBROUTINE MIN4 (PD,AM4,BM4) 

IFIPD. LT .-.425 )G0T010 
IFIPD .GT. . 275 ) GOTO 1 1 
AM 4 =.006 
G0T012 

10 AK4=-.0925*PD-.0335 
GOTO 1 2 

11 AM4 = • 0 7 2 5 * P D- . 0 1 4 

12 IFIPD .GT. .88) G0T013 
IFIPD .LT. -.36) GOTO 14 
BM4=. 184*PD+.076 
RETURN 

13 BM4=.5*PD-.202 
RETURN 

14 BM4= . 4 * PD+ ♦ 145 
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RETURN 

END 

SUBROUTINE MIN8 (PD,AM8,8M8) 
IF(PD .LT. -.55) GO TO 10 
AM8=.07*PD 
GOTO 1 3 

10 AM8=-.028*PD-.054 

13 I F ( PD .LT. - . 45 ) GOTO I I 
I F ( PD .GT. .42) GOTO 1 2 
BM8=.09*PD+.2760 
RETURN 

11 BM8= . 2 I»PD+ . 33 
RETURN 

12 BM8=. 283#PD+. 1652 
RETURN 

END 

SUBROUTINE MI ONE (PD,AMI,BM1) 

AM1=.064*PD-.010 

I F ( PD .GT. . 6) GOTO 11 

IF(PD . LE. .400) GOTO 10 

RM1=-.88*PD+. 743 

RETURN 

10 B M 1 = . 391 
RETURN 

11 B Ml-. 215 



RETURN 

END 
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SUBROUTINE SHAFT ( S N , Q f : ) 

THIS SUBROUTINE COMPUTES SHAFT FRICTION 

ENTERING ARGUMENT IS SHAFT ROTATION RATE (SN) IN REV PER SECOND 
OUTPUT IS FRICTION TORQUE (QF) IN FT-LBS 

S N = A G S ( S N ) 

IFISN.LE. . A 1 7 ) GOTO 1 0 
QF=SN*6000. 

RETURN 
QF = 2 500 . 

RETURN 

END 
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APPENDIX E 
Ship Dynamics 



1. Ship Propulsion Equations 

The ship propulsion equations with no external forces or 
. rudder movement are the two first-order, non-linear, differential 
equations given below: 

ra ~ = ? - R (thrust equation) 



dn 



dt 



Q D’ Q F oQ p 



(torque equation) 



v = ship velocity (ft/ sec) 

T = thrust (lbs) 

R = resistance (lbs) 

2 

ra = mass of ship plus 1G$ (iVrJ££ ) 

i u 

I - rotational inertia of drive shaft (Ib-ft-sec ') 

Qp = prime mover drive torque on propeller (Ib-ft) 

Q„ = shaft friction torque (Ib-ft) 

Qp = propeller torque (Ib-ft) 
n = propeller angular speed (rev/ sec) 

The two equations are coupled through the propeller terns T and 
Qp« The equations were solved simultaneously on an IBM 11 3 G 
co&puter using a fourth order Runge-Kutta integration technique. 

2. Sig n Convention 

The sign convention used here is that propolier angular speed, 
n, and velocity, v, are defined as positive for forward ship motion. 
Torques are considered positive when acting in the direction of 
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APPENDIX E (Cont'd) 



2. (Cont'd) 

positive n; thrusts are considered positive in the direction of 
positive v. 

3- Ship Displacement 

For all simulations, a ship with a displacement of 4000 tons 
was used. In the equations of motion, 10$ of this figure was 
included to account for entrained water. 



APPENDIX F 



FT^A-2 Drive Train 



Digital Simulation of an FT4A~2 Forward Drive Train 

Figure (F»l) is a sketch of a drive train powered by two 

FT4a~ 2 marine gas turbines that was used in the simulations. 

Except for the propeller the drive train is basically the same 

used by Rubis in Reference (11). The reduction gear ratios relate 

propeller speed and torque to power turbine speed and torque. 

The drive train moment of inertia was assumed to be the same 

5 2 

as Rubis*, viz., 2.88 x 10 Ib-ft-sec (referred to propeller 
shaft) . 

It should be noted that at high speeds and low fuel rates, th 
FT4A“2 gas turbine subroutine vjill yield negative torques. A 
statement was made in the main program to ignore this and sub- 
stitute instead a linear windage torque of tho form, 

Q = -2H + 3000 (ft-lb) 

N = Free Turbine Speed (RPM) 

For the crashback simulation, a Falk clutch whose torque 
capacity is 

Q c = 5687.5 P net 

where 

p net = net air pressure acting on shoes 
was used (cf. section 3»2). A supply air pressure of 200 psi was 
assumed. The criteria for clutch disengagement was the same used 
by Rubis in Reference (11): whenever the propeller torque "seen" 

by the clutch exceeded clutch torque, clutch disengagement was 
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1. (Cont f d) 

was assumed to have taken place instantaneously. For the crash- 
back simulation, the equation for f orward clutch disengagement 
torque was 

Q c = 5687.5 (195-30t-4N 2 ) (ft-lb) 
t = time (sec) 

N = propeller shaft rotation rate (rev/sec) 

The propeller torque opposing the forward clutch was 




= reduction gear ratio (see fig. (?•»!)) 

Qp = propeller torque ft-lbs 

During reverse clutch engagement, reverse clutch torque is given 
by 



Q cr = 5687.5 (5t-5-4n 2 ) (ft-lb) 



t = time (sec) 

n = propeller shaft rotation rate (rev/sec) 

The open drive train is defined to exist whenever the quantity 
(5t-5-4n ) is negative; if this occurs, Q„ is set to a zero 

UK 

value obviating a negative value for clutch torque. A slip speed 
across the clutch exists during reverse clutch engagement and 
this is given by 

"slip ■ *3* - 

r v 2 



N - propeller rotation rate (rpra) 
= power turbine speed (rpm) 
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APPENDIX F (Cont’d) 



(Cont*d) 



Kj > K ? , K^, = reduction gear ratios 



The product of N g ^ and is the power dissipated in 
and the integral of power over time is a measure of the 
absorbed in the clutch. 

Reverse clutch engagement was assumed to exist the 
l^slio ec l ua -*- s< * zero. 



the clutch, 
energy 

instant 
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APPENDIX G 

IJ-125Q0-A Drive Train 

1. Digital Simulation of an LH2500-A Powered Drive Train 

Figure (G-l) is a sketch of the drive train powered by two 
LM2500-A tnarino gas turbines with a controllable reversible pitch 
propeller (CRPP). A popular control strategy is to increaso speed 
by first increasing propeller pitch while keeping shaft rotation 
rate constant; once full pitch is attained, shaft rotation rate 
is increased. (2b) A proportional- integral controller of the form 

h f =5 + ^ 

ter, 

Wp = fuel rate (lb/ hr) 

RPM = (ordered shaft RPM)- (Actual Shaft RPM) (RPK) 
s = differential operator 

was used in one simulation to keep shaft rotation rate constant 
while increasing pitch. 

As in the FT4A-2 subroutine, the LM250Q-A subroutine will 
also yield negative torques for high power turbine RPM and low 
fuel rates. These torques were also ignored and a \-dndage torque 
expressed by 

= -1.965 N + 5500 (ft-lb) 

N = power turbine speed (RPM) 
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PAGE 1 

SUBROUTINE EOS IK 

C ACCELERATING FROM 5 KNOTS WITH THE FT4A-2 

COMMON T, OT, Y(20)» DY(20)» STIME, FT I ME » NFWOT 
VA=.96»Y ( 2 ) 

CALL PRPLR (VA,CT,CO) 

DENOM=VA**2+( 33.*Y(1))**2 
C COMPUTE PROPELLOR TORQUE Y(5) 

Y ( 5 ) =262. 5«DENOM«CC 
C COMPUTE PROPELLOR THRUST 

Y ( 6 ) =171 . 5*DENOM*CT 

C PROPELLOR THRUST IS REDUCED BY THRUST DEDUCTION FACTOR 

Y ( 6 ) = Y ( 6 ) * . 9 8 

C ETA IS REDUCTION GEAR EFFICIENCY AS A FUNCTION OF SHAFT SPEED 

E TA=.75+. 22 *( 1.0- EXP ( - Y ( 1 1/1.04 ) ) 

C Y ( 7 > IS FUEL FLOW A A FUNCTION OF TIME FOR THIS MANEUVER 

Y (71=2 500. +( 1 .-EXP ( -T ) ) *9500. 

C RT IS SHIP RESISTANCE AS A FUNCTION OF SHIP'S SPEED 

RT=104.8*(Y(2 ) )«*2 
C SHAFT COMPUTES SHAFT FRICTION 

SN=Y ( 1 ) 

CALL SHAFT (SN,CF) 

AN3=861.*Y( 1) 

WF=Y( 7) 

C FT4A2 COMPUTES POWER TURBINE TORQUE 

CALL FT4A2 (AN3,WF,Q3) 

DELAY=.3+.7*( l.-EXP(-.4*T) ) 

Q3 = Q3-8 14 « 

Q3=Q3»0EL AY 
Q3=Q3+814. 

C Y ( 9 ) IS POWER TURBINE TORQUE 

Y ( 9 ) =Q3 

C QD IS POWER TURBINE TORQUE MODIFIED BY REDUCTION GEAR FACTOR 

0D=28.72*03*E1 A 

C Y ( 1 ) IS SHAFT SPEED IN RPS 

DY ( 1 ) = ( QD-QF-YI 5 ) )»5.54E-7 
C Y ( 8 ) IS SHAFT SPEED IN RPM DIVIDED BY 100. 

Y ( 8 ) = . 6*Y ( 1 ) 

C Y ( 2 ) IS SHIP'S VELOCITY IN FPS 
DY(2)=(Y(6)-RT)^3.33E-6 
C Y ( 3 ) IS DISTANCE TRAVELLED BY SHIP 

DY ( 3 ) =Y ( 2 ) 

C Y ( 4 ) IS SHIP'S VELOCITY IN KNOTS 

Y ( 4 ) = . 59 1 *Y ( 2 ). 

C Y ( 1 6 ) IS PROPELLOR THRUST COEFFICIENT 

Y ( 16) =CT 

C Y ( 1 7 ) IS PROPELLOR TORQUF COEFFICIENT 

Y { 17 ) =CQ 

I F ( T .GT. 10. ) DT = 1 . 

RETURN 

END 



05/03/ 
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'’AGE 1 

SUBROUTINE EQSIK 

C STOPPING THE SHIP WITH THE FTAA-2 

COMMON T, DT, YI20) , DYI20) , S T I ME , F T I ME , NE WOT 
VA=.96*Y( 2) 

CALL PRPLR ( V A t C T i CQ ) 

DEN0M=VA**2+ (33.*Y(1))«*2 
C COMPUTE PKOPELLOK TORQUE Y(5) 

Y(5)=262.5*DEN0K*CQ 
C COMPUTE PROPELLOR THRUST 

Y(6)=171.5*DEN0M*CT 

C PROPELLOR THRUST IS REDUCED BY THRUST DEDUCTION FACTOR 

Y(6)=Y(6)*.98 

C ETA IS REDUCTION GEAR EFFICIENCY AS A FUNTION OF SHAFT SPEED 

ETA=.75 + .22-*( 1 . O-EXP ( - ABS ( Y ( 1) 1/1.04) ) 

C Y ( 7 ) IS FUEL FLOW AS A FUNCTION OF TIME FOR THIS MANFUVER 

Y(7)=16Q0.+26Q0.*( T- 11. 1 19) 

I F ( Y ( 7 ) .GT. 12000. )Y(7)=1?000. 

C RT IS SHIP RESISTANCE AS A FUNTION OF SHIP'S SPEED 

RT=104.8*(Y(2))**2 
C SHAFT COMPUTES SHAFT FRICTION 

$N= Y ( 1 ) 

CALL SHAFT ( S N » Q F ) 

AN3=600.25*(-Y( 1 ) ) 

C Y(9) IS POWER TURBINE SPEFD IN RPM 

Y ( 9 ) = AN3 
W F = Y ( 7) 

C FTAA2 COMPUTES GAS TURBINE TORQUE 

CALL FT4A2 ( AN3, WF, Q3 ) 

I F ( Q3 .LT. 0. )Q3=-2. *AM3+3000. 

Q3=Q3-2165.2 

DELAY=.3+.7*( 1 . -EXP ( - . 4* ( T- 1 1 . 1 19 ) ) ) 

Q3=Q3* DELAY 
03=03+2165.2 

C Y ( 8 ) IS POWER TURBINE TORQUE DIVIDED BY 100 

Y(8)=Q3*.01 

C QD IS GAS TURBINE TORQUE MODIFIED BY REDUCTION GEAR FACTOR 

QD= - 2 0 . 1 *Q3 * E T A 

C Y ( 1 ) IS SHAFT SPEED IN RPS 

DY ( 1 ) = ( QD-OF- Y ( 5 ) )*6.75E-7 
C Y ( 2 ) IS SHIP'S VELOCITY IN FPS 

DY(2)=(Y(6)-RT)*3.33E-6 
C Y ( 3 ) IS DISTANCE TRAVELLED BY SHIP 

DY ( 3 ) = Y (2) 

C Y ( A ) IS SHIP'S VELOCITY IN KNOTS 

Y(4)=.591*Y(2) 

C Y ( 1 9 ) IS THRUST COEFFICIENT 

Y ( 19 ) =CT 

C Y ( 2 0 ) = IS TORQUE COEFFICIENT 

Y ( 2 0 ) = C 0 

I F ( Y ( 2 ) . L E . 0. )FT!ME = T 
I FIT .GT. 16. ) DT= . 4 

RETURN 

END 



05 / 03 ; 
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Y ( 10}=.98*Y( 10) 

RT IS SHIP RESISTANCE IN LBS AS A FUNCTION OF SHIP’S SPEED 
RT= 104, 3*Y(2)*Y(2) 

SHAFT COMPUTES SHAFT FRICTON AS A FUNCTION OF SHAFT SPSEC 
S N = Y ( 1 ) 



144 



5: 

CL 

cl 



o 



CC 

LJ 

U- H- 

3: O 

~ < 

LL O' UL 

O O' 

< a' 

z: *> <1 

o r- uj 

< C3 

[** (*« 

U'O 2 : 

2: < o 

ZD *“ ^ 

U- Uf\ b~ 

< O 

<1 ID 

ST O 

CO <x UJ 

< ^ CC 

03 C 

LU < O >- 

r- o oj 

cj cm in 

o: <1 in o ^ 

0*4- l-U — 

I — vO 2_ **•< v?" 

r-< CL IL o 

< CC e 

UJ ► # O rW 

z co in o s 



a 








M 


c 


O 




X 






r- 


















IXJ 








CO 


< 


O' 






<— ( 




1 


















UJ 








CX 


••* 


• 




UJ 


w > 




UJ 


MJ 


1— 














CL 








D) 


0 


rH 




3) 


> 




h* 


1 


it 


— ^ 






H 






<0 








1 — 


0£ 


I 




O 


1 




m 


LU 


IU 


1 — 






II 
















O 


II 




fX 


w 




• 


ro 


X 


1 






UJ 






UJ 








cx 


H 


O 




0 


CL 




<> 


n 


►— < 








X 






Z 








UJ 


•- 


a: 




h~ 


X 






0 


I'- 


CL 




in 


( — < 






*— 4 








3c 


LL 


O 






UJ 


< 




03 


LL 


X 




« 


h 






CO 








a 


:< 


H- 




UJ 


1 


H 




>r 




UJ 




li 


UL 






0' 








CL 


► 






Z 


0 


UJ 


O' 


< — . 


CO 


•K 




h- 


















X 


t 




KH 


• 




— 


1 — 


« 


» 




0 


• 






K 








CO 


CL 


O 




cn 


f-v 




x 


Ol 


1 


0 






O 














UJ 


CX 






cx 




r—» 


1 






0 




in 








(X 


rH 






H 




• 






•* 


H 


LL 




© 


*-H 






» 






UJ 


*- 






3 




h- 




p 


CM 




CJ 


O 


LU 


C' 




<- 


uu 






Iv 


> 






a 


O 


— J 


0 




(M 


> 


1 


H 




•f 












O 


** 






2 f* 


m 


• 


(V 


(O 


• 




a 


w 


• 


• 


CM 


• 


• 






CL 


0 


2 l 




C* 


cm 




a 


<L 


+ 


CM 


LJ 


>* 




CJ 


• — 


H- 










r— < 


O. 


— * 


CJ 


X 


O 


i— 


0 


IO 


h- 




w 




0 


>' 


O 








CO 


<J 


cn 


vO 






Ct 


11 




r- 


e 


li 


II 


ir 


CO 


II 


« 


CM 


—r 

-U 




K-« 


CO 


n 


w 


0 




0 


<*•* 


CO 


c 


CO 








II 








CL 






1! 




> 


U 


T-l 


H* 


r — < 


< 


1 ! 


CM 


rH 


CM 


> 


- — 


<r 


h~ 


> 


LO 




X 


X 


r- 


II 


cm 


-J 


«— 


• — 1 




<1 


II 


w- 






vO 


w 


— 


w 


h 


CJ 


CL 


O. 


w» 


u. 


z 




LL 


w- 


LJ 


K 


O 


> 


>- 


LL 




>- 


UL 


IL 


LU 


Z 


cl 


(X 


>- 


3: 




O 


» — 1 


>- 


c? 


LU 


O 


O 


O 


M 




0 






DC 


UJ 



o 



o 



o 



APPENDIX H 



KV Major 12 Drive Train 



Figure (H-l) is a sketch of the drive train used for simulation 
2 

The WR of the engine-flywheel was provided by the Hawker-Siddeley Co 
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PAGF 1 

SUBROUTINE EQS1M 

C COASTDOWN AND OPEN CLUTCHING WITH THE MIRRLESS ENGINE 

COMMON T, DT , Y ( 20 ) , DY (20) , STIME, FTIME, NEWDT 

Y ( 1 ) IS SHAFT SPEED IN RPS 
Y ( 2 > IS SHIP'S SPEED IN FPS 
Y ( 3 ) IS DISTANCE TRAVELLED IN FT 

Y ( A ) is PROPELL.OR THRUST 
Y ( 5 ) IS PROPELLOR TORQUE 

Y A= . 96* Y ( 2 ) 

CALL PRPLR (VA,CT,CQ) 

DENOM=VA**2 + ( 33 .*Y ( 1 > > **2 
COMPUTE PROPELLOR TORQUE 
Y ( 5 ) = 262 . 5*DENOM*CQ 
COMPUTE PROP C LL OR THRUST 
Y(A)=171.5*DE NOM * C T 

PROPELLOR THRUST IS REDUCED BY THRUST DEDUCTION FACTOR 
Y ( A ) = .98*Y ( A ) 

ETA IS REDUCTION GEAR EFFICIENCY AS A FUNCTION OF SHAFT SPEED 
ETA=.75 + .22* ( 1 .O-FXP (-Y(l) / 1 .OA ) > 

Y 6) IS SHIP'S SPEED IN KNOTS 
Y ( 6 ) = . 59 1 * Y ( 2 ) 

Y ( 7 ) IS FUEL SETTING IN PER CENT 
Y(7)=97.B-20.*T 
I F ( Y ( 7 ) .LE. I2.3)Y(7) = I2.3 

TORQ COMPUTES MIRRLESS DIESEL TORQUE AS A FUNCTION OF 
FUEL SETTING 
T 0RQ=550 . * Y ( 7 ) 

C RT IS SHIP RESISTANCE AS A FUNCTION OF SHIP'S SPEED 

RT= 1 OA . 8 * Y ( 2) #*2 

C SHAFT COMPUTES SHAFT FRICTION 

S N = Y ( 1 ) 

CALL SHAFT (SN,OF) 

C Y ( 8 ) IS DIESEL TORQUE 

Y ( 8 ) =TORQ 

C QD IS DIESEL TORQ MODIFIED BY REDUCTION GEAR 

QD=3.72*TORQ*ETA 
DY ( 1 ) = ( QD-Y ( 5 ) ) *-8 . 6AE-6 
DY(2)=(Y(4)-RT)*3.33E-6 
DY ( 3 ) = Y ( 2 ) 

C QCF IS FORWARD CLUTCH TORQUE 

QCF=60000.*EXP(-T/5.0) 

C QSLK IS PROPELLOR TORQUE AS SEEN DY CLUTCH 

OSLK= Y ( 5 1/3.72 

I F ( APS ( QCF ) .LT. ABS ( QSLK ) > FT IME = T 

Y ( 9 ) =QC F 

Y( 1 0 ) =QS L K 

RETURN 

END 
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SUBROUTINE EQSIM 
KIRRLFSS ENGINE 

OPEN-CLUTCH — RECLUTCHING PHASE 

COMMON T, DT, Y ( 20 ) , DY(20) , ST I ME , F T I ME , NE WDT 

Y ( 1 ) IS SHAFT SPEFD IN RPS 

Y ( 2 ) IS SHIP'S VELOCITY IN FPS 

Y ( 3 > IS DISTANCE TRAVELLED BY SHIP IN FEET 

Y ( 4 ) IS ENGINE SPEED IN RPS 

Y ( 5 ) IS CLUTCH ABSORBED ENERGY 

Y ( 6 ) IS PROPELLOR THRUST 

Y ( 7 ) IS FUEL SETTING IN PER CENT 

Y ( 8 ) IS SHAFT SPEED IN RPH 

Y ( 9 ) IS DIESEL TORQUE 

Y( 10) IS CLUTCH SLIP SPEED 

Y ( 1 1 ) IS REVERSE CLUTCH TORQUE 

Y ( 1 2 ) IS POWER ABSORBED BY REVERSE CLUTCH 

Y ( 1 3 ) IS PROPELLOR TORQUE 

V A=.96*Y( 2 ) 

CALL PRPLR! VA,CT,CQ) 

DENON = V A* VA+ ( 33.*Y( 1 ) )**2 
Y ( 13) =262.5*DFN0M*CC 
Y(6)=171.5*DEN0M*CT 

RT IS SHIP RESISTANCE AS A FUNCTION OF SHIP'S SPEED 
RT=104.8*Y( 2) *Y(2) 

SHAFT COMPUTES SHAFT FRICTION 
SN= Y ( 1 ) 

CALL SHAFT! SN,QF) 

R P M = Y ( 4 ) *60. 

Y( 7) =97.8-20. *(T-11. 199) 

I F ( Y ( 7 ) .LE. 12.35Y! 7 ) = 1 2 . 3 
T0RQ=550 . *Y ( 7 ) 

Y ( 9 ) = TORQ 

TORQ COMPUTES MIRRLESS DIESEL TORQUE AS A FUNCTION 
OF FUEL SETTING 

ETA IS REDUCTION GEAR EFFICIENCY AS A FUNCTION OF SHAFT SPEED 
ETA=.75+.22* ( 1.0- EXP ( - Y ( 1 ) /I .04 ) ) 

CLCRT IS CLUTCH REVERSE TORQUE 
CLCRT = 60000 . * ( 1 . - EXP ( ( - T + 1 1 . 1 99 ) / 5 . ) ) 

DY ( 1 )=-CLCRT*6.8!5E-5-.921E-5*QF-1.842E-5*Y( 13) 

DY(2)=( Y(6)-RT)*3.33E-6 
DY(3)=Y(2) 

DY(4)=(T0RQ-CLCRT)*2.21E-4 
Y( 1 0 ) — ( 3.72*Y( l)-(-RPM) >*60. 

Y( 1 2 ) = Y ( 10) *CLCRT*1 .904E-4 
DY(5)=5'50.*Y( 12) 

Y ( B ) =60 . *Y ( 1 ) 

Y ( 1 4 ) IS SHIP'S SPEED IN KNOTS 
Y(14)=.591*Y(2) 

I F ( Y ( 10) .LE. 0. ) F T I M E = T 

RETURN 

END 
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SUBROUTINE EG'SIM 

C STOPPING THE SHIP WITH KIRRLESS 

COMMON T,DT,Y(20) »DY(20) , ST I ME, FT I ME, NFWDT 
Y( I ) IS SHAFT SPEED IN RPS 
Y ( 2 ) IS SHIP'S VELOCITY IN FPS 
Y { 3 ) IS DISTANCE TRAVELLED IN FT 
Y(4) IS PROPELLOR TORQUE 
Y { 5 ) IS PROPELLOR THRUST 
Y { 6 ) IS FUEL SETTING IN PER CENT 
VA=.96*Y ( 2 ) 

CALL PRPLR (VA,CT,CQ) 

DENOM = VA*VA + ( 33-*Y( L ) )«*2 
Y { A )=262. 5*DENOM*CQ 
Y(5)=171.5*DEN0M*CT 
RT IS SHIP RESISTANCE 
RT=104.8*Y( 2 ) *Y( 2 ) 

SHAFT COMPUTES SHAFT FRICTION 
S N = Y ( 1 ) 

CALL SHAFT ( S N , Q F ) 

ETA IS REDUCTION GEAR EFFICIENCY AS A FUNCTION OF SHAFT SPEED 
ETA=.75+.22*{ I . O-EXP ( - ABS ( Y ( I ) 1/1.04) ) 

Y{6)=12.3+15.*(T-17. 199) 

I F { Y { 6 ) .GF. 97.6 )Y( 6 1=97.8 

TORO COMPUTES MIRRLESS DIESEL TORQUE AS A FUNCTION 
OF FUEL SETTING 
T0RQ=550.*Y ( 6 ) 

Y ( 7 ) = TORQ 

DELAY=I.-5XP{-(T-I7. 199)/. 15) 

Y(7)=Y(7)-6765. 

Y{ 7 )=Y{ 7 ) »DFL AY 
Y ( 7 ) =Y ( 7 ) +6765. 

C QD IS DIESEL TORQUE MODIFIED BY REDUCTION GEAR 

GD=-3.72*T0RQ*BTA 
C Y ( 8 ) IS DIESEL RPM 

Y(8)=223.*(-Y( 1) ) 

C Y ( 9 ) IS BEARING LOAD TORQUE LIMIT 

Y(9)=20400.+88.6*Y(8 ) 

DY ( 1 ) = (0D-Y(4) )*8.64E-6 

C THRUST IS REDUCED BY THRUST DEDUCTION FACTOR 

Y(5) = .98*Y< 5) 

DY ( 2 ) = ( Y { 5 ) -RT ) «- 3 . 33 E-6 
DY ( 3 ) =Y ( 2 ) 

C Y ( 10)' IS SHIP'S SPEED IN KNOTS 

Y ( 10)=.591*Y( 2) 

I F ( T .GT. 20. )DT=1. 

I F ( Y ( 2 ) . L E . 0 . ) F T I M F = T 

RETURN 

END 



05/03/ 



150 



APPENDIX I 

D & W 7K98FF Drive Train 

Figure (I-l) is a sketch of the 7K98FF drive train. The WR 2 
of the engine v?as not available in the literature surveyed. An 
official of the Burmeister and Wain Corporation provided the following 
weights : 

Piston ..... 4.18 tons 
Crosshead... 3*35 tons 
Crank 4.97 tons 

From the "Shock and Vibration Handbook", Vol. 3, p. 38~2» the following 

p 

expression was found for finding the WR' of a diesel: 



• / W 

J = JL 



+ W Q (1- |- ) ] R 



W p = weight of piston, piston pin, and cooling fluid (lb) 

W - weight of connecting rod (lb) 

h = fraction of rod length frora the lower center of gravity 
R = crank radius (in) 

From the above expression, the polar moment of inertia was estimated 
at 4.25 x 10 3 lb-ft-sec 2 . 
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